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ABSTRACT 

The  simulation  program  GCAAS  (General  Circular  Antenna  Array  Simu¬ 
lator)  has  been  developed  to  simulate  the  output  data  from  a  circularly 
disposed  antenna  array  (CDAA)-direction  finding  system  under  known  input 
signal  conditions.  The  CDAA  system  which  has  been  simulated  is  similar 
to  the  Wullenweber  Direction  Finding  System  at  the  University  of  Illinois 
and  may  consist  of  as  many  as  360  antenna  elements  in  a  circular  array  of 
any  reasonable  dimensions  and  simulates  the  outputs  obtained  from  the  array 
used  in  a  scanning  mode. 

The  input  signal  conditions  may  include  up  to  five  incoming  signals 
at  different  azimuthal  and  elevation  angles  of  arrival  with  time-dependent 
relative  phase  shifts  between  the  incoming  signals.  The  output  consists 
of  the  sum,  difference  and  differential  phase  data  sampled  at  regular  inter¬ 
vals  of  azimuth  and  stored  on  magnetic  tape  for  further  use  in  the  deter¬ 
mination  of  the  direction  of  arrival.  A  maximum  of  100  samples  of  the 
output  data  set  may  be  obtained  from  each  scan. 
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I.  INTRODUCTION 

This  generalized  simulator  of  a  circularly  disposed  antenna  array 
(CDAA)  direction  finding  system  has  been  developed  to  satisfy  the  need  for 
a  means  of  producing  simulated  system  output  data,  under  known  input  con¬ 
ditions,  to  be  used  in  future  studies  relating  to  the  different  methods  for 
direction  of  arrival  calculation.  As  a  consequence  of  the  general  nature  of 
the  simulator,  it  may  also  be  used  in  studies  related  to  the  construction 
and  operation  of  the  CDAA  system. 

The  system  outputs  simulated  are  the  scanning  mode  sum,  difference* and 
differential  phase  information;  it  is  possible  in  this  simulator  to  produce 
these  outputs  corresponding  to  a  series  of  related  time-dependent  scans. 

The  primary  needs  for  thi9  simulated  data  are  for  u9e  in:  1)  development  of 
a  difference  method  for  direction  of  arrival  calculation,  2)  a  study  of  the 
presently  used  sum  and  differential  phase  methods  for  direction  of  arrival 
calculation,  and  3)  a  comparison  of  the  three  methods  for  direction  of  arrival 
calculation. 

The  known  input  conditions  this  simulator  has  been  designed  to  handle  are: 

1)  Single  input  signal. 

2)  Multiple  input  signals  from  the  same  source  (wave  interference). 

3)  Amplitude  modulation  of  the  signal. 

4)  Addition  of  random  noise  to  the  signal. 

The  individual  signal  conditions  which  comprise  the  total  input  signal 
conditions  are  to  be  specified  (for  each  input  signal)  by  azimuthal  and 
elevation  angle,  relative  amplitude,  initial  reference  phase  and  a  time- 
dependent  reference  phase  shift. 
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II.  THE  PHYSICAL  SYSTEM 

The  CDAA  direction  finding  system  to  be  simulate!  is  similar  to 
the  Wullenweber  Direction  Finding  System  at  the  University  of  Illinois. 

The  array  consists  of  a  number  of  identical  folded  vertical  monopole 
antenna  elements,  mounted  on  a  ground  mat,  evenly  spaced  around  the  cir¬ 
cumference  of  a  circle.  Within  this  circular  array  and  concentric  to  it 
is  located  a  vertical  reflecting  cylinder  of  a  height  somewhat  greater 
than  that  of  the  antenna  elements.  The  output  of  each  antenna  element 
is  connected  by  a  matched  cable  to  a  multicoupling  device;  each  cable 
has  the  9ame  electrical  length.  Each  multicoupler  iB  a  wide  band  ampli¬ 
fier  with  a  gain  of  unity  and  several  outputs  with  the  same  output  signal 
available  at  each  of  the  multicoupler  outputs.  (Refer  to  Figure  1.) 

One  output  from  each  multicoupler  ,1s  fed  into  an  electro-mechanical 
scanner.  Thig  scanner  selects  a  group  of  consecutive  antenna  elements  to 
produce  a  number  of  scanning  element  voltages.  These  scanning  elements 
have  the  same  angular  spacing  as  the  antenna  elements  and  are  subdivided 
into  a  right  and  left  bank,  each  of  which  has  the  same  number  of  elements. 
Midway  between  the  two  banks  is  a  plane  of  symmetry  referred  to  as  the 
boresight.  The  scanner  is  also  set  to  rotate  at  a  constant  speed  in  the 
clockwise  direction.  Thus,  the  scanner  is  continually  selecting  new 
anta-ra  elements  to  produce  the  scanning  element  voltages  as  the  boresight 
moves  to  the  right.  (Note:  The  clockwise  direction  is  referred  to  as 
the  right  and  the  counter-clockwise  direction  as  the  left.) 

The  right  and  left  scanning  elements  are  connected  to  the  inputs  of 
two  identical  sets  of  delay  lines,  one  set  for  each  bank.  The  effect  of 
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ANTENNA  ELEMENTS, 
CABLES  8  MULTICOUPLERS 


Figure  l.  Flock  Diagram  of  CDAA  Direction  Finding  system. 
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Che  addition  of  the  delav  lines  is  to  make  the  delayed  scanninp,  element 
voltages,  for  an  incoming  signal  arriving  at  a  particular  elevation 
angle  of  arrival  known  as  the  cophasal  angle,  appear  to  be  the  same  as 
the  element  voltages  that  would  be  received  from  a  linear  array  of  verti¬ 
cal  antenna  elements  rotating  about  a  vertical  axis  at  the  center  of  the 
arrav . 

The  outputs  of  the  delay  lines  are  tied  to  identical  summing  net¬ 
works,  one  for  each  bank.  These  summing  networks  perform  the  phasor 
addition  of  the  delayed  scanning  element  voltages  that  are  the  inputs  to 
the  network.  The  outputs  of  the  summing  networks  are  the  right  and  left 
bank  voltages. 

It  is  from  these  two  bank  voltages  that  the  system  outputs  are  ob¬ 
tained.  Feeding  the  bank  voltages  into  the  hybrid  network,  the  sum  and 
difference  voltages  are  obtained;  feeding  the  bank  voltages  into  the 
phase  meter  yields  the  differential  phase  data. 

The  system  produces  continuous  analog  data, but  for  acquisition  and 
storage  of  the  data  in  a  digital  format  the  data  must  be  in  a  discretely 
sampled  format.  The  data  acquisition  system1  performs  the  necessary  sam¬ 
pling  and  analog- to-digital  conversion  to  put  the  output  In  a  format  suit¬ 
able  for  tape  storage. 

The  data  acquisition  system  Is  activated  by  the  position  of  the 
scanner  boresight.  Once  the  acquisition  system  is  activated,  the  sum, 
difference, and  differential  phase  data  are  sampled  and  converted  from 
voltage  level  to  digital  representations  and  stored  on  magnetic  tape. 

The  activation  settings  in  the  scanner  are  such  that  a  given  number  of 
sample  points  per  scan  is  taken  with  equal  angular  intervals  between 
the  boresight  positions  for  consecutive  sample  points.  The  entire  angu¬ 
lar  interval  over  which  the  sample  points  are  taken  is  referred  to  as  the 
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sampling  window.  The  window  Is  always  less  than  90°  and  the  center  of 
the  window  has  been  designated  the  sector  center. 

It  is  from  a  later  analysis  of  the  Information  stored  on  the  magnetic 
tape  that  the  direction  of  arrival  Is  calculated. 
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III.  MATHEMATICAL  MODEL  OF  THE  SYSTEM 

The  model  of  the  CDAA  consists  of  NOA  antenna  elements  located  out¬ 
side  a  reflecting  cylinder  of  radius  A  with  the  distance  from  the  antenna 
element  to  the  reflecting  cylinder  surface  equal  to  D. 

In  the  model,  an  antenna  element,  cable,  and  multicoupler  are  con¬ 
sidered  one  unit.  Since  all  or  the  cables  are  of  the  same  electrical 
length  and  the  gain  of  each  multicoupler  is  unity,  the  outputs  of  the 
multicouplers  are  taken  to  be  the  same  as  the  outputs  of  the  antenna 
elements  to  which  they  are  connected.  The  phasor  voltage  output  cal¬ 
culations  used  for  the  vertical  folded  monopole  antenna  elements, 
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excited  by  a  single  signal  are  basically  the  same  as  those  used  by  Jones. 
These  calculations  have  been  expanded  for  multiple  signals  (wave  inter¬ 
ference)  by  the  addition  of  a  relative  amplitude  factor  and  a  reference 
phase  for  each  signal. 

Jones  '  modified  expressions  for  antenna  element  voltage,  as  modified 
for  element  i  and  signal  w,  are: 

Amplitude 

A.  eJ  *  [c°S(e  )  *  SIN  2tt  *  D  ,,  C0S(O  *  COS(e)  }  *  WAMP  (1) 

1W  1W  W  W  Vrr  1W  W  W 


Phase 


Jiw  (<t>iw’ 


0  ) 
w 


2tt  *  A 
WL 


*  cos  (<p.  )  *  cos  (e  )  +  wphi 

*  iw  w  w 


(2) 


where 


WL  *»  Wavelength  of  the  signal,  calculated  from  the  carrier  frequency 
FREQC  (in  MHz)  as  WL  -  299.8/FREQC 
A  ■  Radius  of  the  reflecting  cylinder  in  meters 
D  ■  Antenna  element  to  reflecting  cylinder  distance  in  meters 
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a  =  Difference  between  the  azimuthal  angle  of  the  antenna  element 

iw 

i  and  the  angle  of  arrival  of  signal  w 

0  -  Elevation  angle  of  signal  w 

w 

WAMP  =  Signal  w  relative  amplitude  factor 
w 

WP111  =  Signal  w  reference  phase  factor* 

w 

Figure  2  shows  the  geometry  that  relates  these  terms  and  Equations  (1)  and  (2; 

For  the  multiple  input  signal  conditions  the  antenna  element  voltage, 
by  superposition,  is  taken  as  the  sum  of  the  voltages  excited  by  the 
individual  signals.  Since  summations  are  more  readily  carried  out  in 
Cartesian  coordinates  rather  than  polar  coordinates,  the  phasor  voltages, 
which  are  excited  by  the  individual  signals,  are  converted  from  polar 
to  Cartesian  form  and  summed  over  the  individual  incoming  signals  to 
give  the  antenna  element  voltage.  Thus  for  n  incoming  signals  the  com¬ 
plex  antenna  element  voltage,  for  element  i,  is  computed  as: 


REi 

*  W=1 

[A, 

iw 

(*iw’ 

e  )  *  cos 

w 

<pt„ 

(*1W  V>] 

(3) 

IMi 

n 

"  Z  . 

W=1 

[A. 

iw 

(<l>iw’ 

0  )  *  SIN 
w 

(Piv 

ew»>- 

(O 

The  superposition  principle  is  also  used  to  add  the  effect  of  the 
time- dependent  phase  shift,  for  the  individual  signals,  to  the  antenna 
element  voltage.  The  antenna  element  voltage  is  recalculated  for  the 
Time  T  by  adding  a  phase  shift  voltage,  for  each  incoming  signal,  to 
the  element  voltage  at  time  To.  Thus  for  element  i  at  time  T,  and  n 
incoming  signals,  the  element  voltage  is  recalculated  as: 

REi(T)  *  ^i (To)  +  fAlw  (*i„’  V  *  (C0S  (PA2)-C0S(PA1))]  (5) 
IMi(T)  "  IMi(To)  +  tAiw  (^iw’  V  *  (SIN  <PA2>-SIN  (PA1))], (6) 


where 
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PAl  -  P  (ii  w)  +  Time  Dep  phase  shut  at  To,  lor  signal  w  (7) 
lw  lw 

PA2  -  P  U  ,  w)  +  Time  Dep  phase  shift  at  T,  for  signal  w.  (8) 
lw  lw 

For  the  calculation  of  the  time  dependent  phase  9hlft  and  ’he 
amplitude  modulation  of  the  incoming  signals, it  is  necessary  to  compute 
simulated  time.  This  simulated  time  has  as  its  basis  the  rotation  of 
the  scanner.  The  scanner  rotates  once  every  SCANT  seconds  and  if  DEG 
is  the  angle  of  rotation  of  the  boresight  since  T*0,  then  the  simulated 
time  T  may  be  given  by: 

T  -  SCANT  *  (DEG/360)  (9) 

Where  PST^  is  the  amount  of  phase  shift  per  scan  for  signal  w, 
Equations  (7)  and  (8)  become 

PAl  =  P,  («(>.  ,  0  )  +  PST  *  (To/SCANT)  (10) 

iw  Tiw’  w  w 

PA2  -  Piw  (*iw,  0w)  +  PSTw  *  (T/SC/iNT)  •  (11) 

The  scanner  consists  of  NAS  b«.  .ing  elements  with  NAS2  »  NAS/2 
elements  in  each  bank.  The  scanner  will  be  set  to  take  ISAM  sample 
points  per  scan;  the  interval  between  the  boresight  for  consecutive 
sample  points  in  the  scan  is  SAMI.  The  sector  center  is  set  at  SECT 
and  the  width  of  the  sampling  window  will  be  WIDTH  ■  (ISAM-1)  *  SAMI. 

The  scanner  reactively  couples  the  scanning  elements  to  the  antenna 
elements  to  produce  the  scanning  element  voltage.  The  scanner  is  of 
such  a  nature  that  its  coupling  can  be  simulated  by  a  direct  (maximum 
to  zero)  relationship  with  the  displacement  of  the  boresight.  When  the 
scanning  elements  on  the  rotor  are  aligned  with  the  antenna  elements  on 
the  stator  (or  in  other  words,  the  boresight  is  exactly  midway  between 
two  antenna  elements),  the  scanning  elements  will  have  the  same  voltages 
as  the  antenna  elements  with  which  they  are  aligned  and  the  adjacent 


•  • 


10 

antenna  elements  will  have  no  effect  upon  the  scanning  element  voltages. 

As  the  rotor  elements  move  out  of  direct  alignment  with  the  stator 
elements,  the  scanning  element  voltages  will  be  taken  as  the  sum  of  the 
effects  of  the  two  adjacent  stator  element  voltages  (See  Figure  3.) 
Consider  a  stator  element  1  which  is  to  the  left  or  aligned  with  the 
rotor  element  i  and  the  stator  element  r  which  is  the  first  stator  element 
to  the  right  of  element  1.  The  displacement  between  adjacenr  stator 
elements  i9  ANOA  *  (360/NOA)  and  if  the  displacement  between  stator 
element  1  and  rotor  element  i  is  CDEL,  the  scanning  element  voltage  will 
be  taken  as : 


SRE  -  [ARE  *  (AN0A-CDEL)/AN0a] 

1  1  +[AREr  *  CDEL/ANOA] 

(12) 

SIM  -  [  AIM  *  (ANOA-CDEL) /ANOA] 

1  1  +[AIMr  *  CDEL /ANOA] 

(13) 

Since  the  scanner  coupling  is  linearly  related  to  displacement, 
it  is  possible  to  have  sampling  intervals  which  may  have  values  other 
than  multiples  of  an  antenna  spacing  interval. 

In  modeling  the  system,  the  addition  of  signal  amplitude  modulation 
and  random  noise  signals  is  most  effectively  .done  at  the  time  the 
scanning  element  voltages  are  calculated.  The  amplitude  modulation 
is  added  by  multipling  the  antenna  element  voltages  by  an  amplitude 
factor  calculated  for  the  sample  point  of  interest  at  that  time,  where 
it  is  considered  that  the  sampling  time  is  relatively  short  in  com¬ 
parison  to  the  modulation  period.  The  percent  modulation  is  defined  as 
PCAM  and  the  modulating  frequency  is  FREQM  (in  KHz).  Using  the  simulated 
time  T,  for  the  sample  point  of  interest,  the  amplitude  modulation  fac¬ 


tor  is: 


Figure  3.  Scanner  Geometry. 


AMS  -  (1-PCAM)  f  PCAM  *  C0S(2r  *  FRF.OM  *  1000  *  T)  .  (14) 

The  noise  signals  are  produced  by  the  use  of  a  random  number  generator. 
The  random  noise  voltage  (ANRE,  ANIM)  added  to  the  antenna  element 
voltages  when  the  scanning  element  voltage  is  calculated,  Is  determined 
by  multipling  the  relative  noise  level  amplitude  ANOISJ  by  successive 
real  numbers  supplied  by  a  pseudo-  random  number  generator.  With  the 
addition  of  amplitude  modulation  and  random  noise  signals,  Equations 
(12)  and  (13)  become 

SREi  -  [(AMS  *  ARE1  +  ANRE^)  *  (ANOA-CDEL) /ANOA] 

+  [  (AMS  *  ARE^  +  ANREr)  *  CDEL/ANOA]  (15) 

SIM  -  [(AMS  *  AIMX  +  ANIM1)  *  (ANOA-CDEL) /ANOA] 

■+■[  (AMS  *  AIM  +  ANIM  )  *  CDEL/ANOA].  (16) 

r  r 

The  delays  which  should  be  introduced  by  the  phase  delay  lines  are 
calculated  from  the  antenna  array  geometry;  for  these  calculations  it  is 
assumed  that  the  scanning  elements  are  aligned  with  the  antenna  elements 
around  the  circumference  of  the  array  The  phase  delay  that  will  be 
introduced  for  element  i,  which  is  i  elements  from  the  boresight,  is 
directly  proportional  to  the  distance  of  signal  propagation,  along  the 
boresight,  from  element  I  to  the  furthest  element  from  the  boresight, 
where  the  signal  azimuthal  angle  of  arrival  is  taken  as  the  boresight 
setting.  This  distance  can  be  realized  as: 

A  *  COS  (CPHI)  *  (C0S((i-l)  *  ANOA  +  AN0/.2) 

-COS  ( (NAS2-1)  *  ANOA  +  AN0A2) ) 

where  AN0A2  •  AN0A/2  and  CPHI  is  the  cophasal  or  assumed  elevation  angle 

of  arrival. (See  Figure  4.)  The  phase  delay  is  then  determined  in  radians 

as:  PDELAY  -  (2-  *  A/WL)  *  COS  (CPHI)  *  (C0S((i-l)  *  ANOA  +  AN0A2) 

-COS ( (NAS2-1)  *  ANOA  +  AN0A2))  .  (17) 
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In  modeling  the  system  it  is  possible  to  combine  the  addition  of 
the  phase  delays  to  the  scanning  element  voltages  and  the  operation  of 
the  summing  network  so  chat  the  bank  voltages  may  be  calculated  directly 
from  the  scanning  element  voltages.  The  phase  delays  are  added  to  the 
scanning  element  complex  voltages,  to  obtain  the  delayed  scanning  ele¬ 
ment  complex  voltages,  as  if  the  phase  delays  were  a  rotation  of  the 
axes  in  a  Cartesian  coordinate  system.  The  summing  network,  as  the  name 
implies,  produces  the  complex  voltage  for  the  bank  by  computing  the 
complex  sum  of  the  delayed  scanning  element  voltages.  Thus  the  combined 
effect  of  the  delay  lines  and  the  operation  of  the  summing  network  will 

determine  the  bank  voltage  as: 

NAS  2 

SUMI  =  Z  [SIM,  *  COS  (PDELAY  )  -  SRL  *  SIN  (PDELAY . ) ]  (18) 

i-1  1  11  1 

NAS  2 

SUMR  =  I  [SIM,  *  SIN  (PDELAY.)  +SRE  *  COS  (PDELAY.)]  .  (19) 

i=l  1  11  i 

The  system  outputs  are  obtained  as  the  outputs  of  the  hybrid  net¬ 
work  and  the  phase  meter  with  the  right  and  left  bank  voltages  used  as 
the  inputs. 

The  outputs  of  the  hybrid  network  are  the  complex  sum  and  difference 
of  the  bank  voltages;  the  magnitudes  of  these  two  complex  quantities  are 
the  sum  and  difference  data  outputs  for  the  system.  These  two  operations 
are  carried  out  as : 

DATAST  =  SQRT  ((SUMRr  +  SUMR^  **  2+(SUMIr  +  SUMl^  **  2)  (20) 

for  the  sum  data,  and 

DATADT  -  SQRT  ((SU1®,  -  STjMR^  **  2+(SUMIr  -  SUMI^  **  2)  (21) 

for  the  difference  data. 

The  differential  phase  measured  by  the  phase  meter  provides  the 
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differential  phase  data  output  and  Is  computed  as 

DATP  =  ATAN2  (SUMI^  SUMR  )  -  ATAN2(SUMI1,  SUMR^  (22) 

Where  ATAN2  is  tan  ^  (SUMI/SUMR)  defined  from  -180  to  +180  degrees. 

In  modeling  the  system  on  a  digital  computer  the  function  of  the 
data  acquisition  system  has  been  simulated,  as  the  output  is  already 
in  a  digital  form  and  can  be  stored  on  magnetic  tape  in  a  digital  for¬ 
mat.  It  should  be  noted  that  in  modeling  the  system  the  relationship 
between  the  position  of  the  scanner  boresight  and  the  simulated  time, 
for  a  given  sequence  of  scans  or  sample  points,  is  the  same  as  the  re¬ 
lationship  between  the  position  of  boresight  of  the  physical  system  and 
real  time,  for  the  same  sequence. 
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[V.  THE  SIMULATION  PROGRAM  -  OPERATION  OF  SIMULATION 

The  simulation  program  has  been  named  GCAAS,  standing  for  General 
£lrcular  Antenna  Array  Simulator,  a  listing  oi  this  program  is  contained 
in  Appendix  A  of  this  report 

In  an  attempt  to  facilitate  a  rapid  understanding  of  the  basic 
operation  of  the  Simulation  Program,  the  program  has  been  reduced,  in 
Reduced  Program  listing  number  one,  to  the  basic  control  statements  and 
input  operations  and,  in  Reduced  Program  listings  two  thru  eight,  to 
the  FORTRAN  statements  used  in  realizing  the  mathematical  model  of  the 
system.  These  reduced  program  listings  are  contained  in  Appendix  B  of 
this  report.  All  of  the  other  program  operations,  which  are  not  in¬ 
volved  in  simulation  calculations,  have  been  represented  in  the  reduced 
listings  as  blocks. 

Before  explaining  the  operations  of  the  program,  it  will  be  helpful 
to  define  certain  variables  and  constants  used  in  the  program. 

Reference  Numbers 

SETNUM  is  the  number  used  to  identify  a  given  data  set  on  the  output 
data  tape,  where  a  data  set  is  the  output  data  stored  on  the  tape  for 
one  execution  of  the  program. 

I RUN  is  the  number  used  to  identify  a  given  data  run,  which  is  one  of 
the  NRUN  data  runs  in  the  data  set  and  where  a  data  run  is  one  simulation 
of  a  given  antenna  array  under  specified  input  and  system  operation 
conditions . 

INAVE  is  the  number  used  to  identify  a  given  incoming  signal  condition, 
where  there  are  NWAVE  incoming  signals  in  the  data  run. 

is  the  number  used  to  Identify  a  given  scan,  where  there  are  ISCAN 


16 


scans  in  the  data  tun 

IP  13  the  number  used  to  identity  a  given  sample  point,  where  there  are 
ISAM  sample  points  in  each  scan  of  the  data  run 

1ANT  Is  the  antenna  element  reference  number,  where  element  number  one 
is  defined  as  being  at  an  azimuth  of  2ero  degrees  and  the  reference 
numbers  increase  in  the  clockwise  direction  up  to  element  number  NOA, 
which  is  the  first  element  to  the  left  of  element  number  one 
I AW  is  the  number  of  elements  an  antenna  element  is  located  away  from 
a  given  incoming  signal  azimuthal  angle  oi  arrival.  (See  Figure  A). 

ISE  is  the  scanning  element  reference  number,  where  element  number  one 
is  the  element  closest  to  the  boresight  and  element  number  NaS2  is  the 
furthest  element  from  the  boresight. 

IDIR  is  used  to  specify  if  1AW  or  ISE  is  to  the  left  (=>1)  or  the  right 
(=2)  of  the  azimuthal  angle  of  arrival  cr  the  boresight,  respectively . 

It  i9  also  used  to  distinguish  between  the  left  (=1)  and  right  (=2) 
scanning  banks, 

K  is  the  reference  number  of  the  antenna  element  to  the  left  of  a  given 
incoming  signal  azimuthal  angle  of  arrival. 

KC  is  the  reference  number  of  the  antenna  element  to  the  left  of  the 
first  scanning  element  in  the  left  scanning  bank  for  a  given  sample 
point. 

IL  is  the  reference  number  of  the  antenna  element  to  the  left  of  scanning 
element  ISE  for  a  given  sample  point 

IR  is  the  reference  number  of  the  antenna  element  to  the  right  of  scanning 
element  ISE  for  a  given  sample  point 
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Figure  4.  Reference  Numbers. 
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Natural  Constants 

VEC  is  the  velocity  of  propagation  (=299.8  meters/microsecond)  . 

PI  is  set  equal  to  the  value  n  (=3.14159265). 

TW0P1  is  9et  equal  to  the  value  2n. 

CONRAD  is  the  radian  to  degrees  conversion  factor  (=360/2it). 

Input  Variables 

ITAPE  is  the  tape  drive  number  (=1,  2,  3)  on  which  the  output  data  tape 
is  mounted. 

MODE  specifies  if  a  new  (=1)  or  c.ld  (=2)  output  data  tape  is  being  used. 
NRUN  is  the  number  of  data  runs  in  the  data  set. 

NOA  is  the  number  of  antenna  elements  in  the  array  for  a  given  data  run. 
A  is  the  reflecting  cylinder  radius  (in  meters)  for  a  given  data  run. 

I)  is  the  reflecting  cylinder  to  antenna  element  spacing  (in  meters)  for 
a  given  data  run. 

NAS  is  the  number  of  scanning  elements  for  a  given  data  run. 

CPHI  is  the  cophasal  angle  (in  degrees)  for  a  given  data  run. 

NWAVE  is  the  number  of  incoming  signals  for  a  given  data  run. 

FREQC  is  the  carrier  frequency  of  the  incoming  signals  (in  MHz)  for  a 
given  data  run. 

FREQM  is  the  amplitude  modulation  frequency  of  the  incoming  signals  (in 
KHz)  for  a  given  data  run. 

PCAM  is  the  amplitude  modulation  factor  for  the  incoming  signals 
(0  £■  PCaM  s:  1.0)  for  a  given  data  run. 

WAMP  is  the  relative  RMS  amplitude  for  a  given  incoming  signal  condition. 
WPHI  is  the  reference  phase  (in  degrees)  for  a  given  incoming  signal 
condition. 

AZM  is  the  azimuthal  angle  of  arrival  of  a  given  incoming  signal 
(0  <  AZM  <  360  degrees). 
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ELF.V  is  the  elevation  angle  of  arrival  lor  a  given  incoming  signal 
(0  *  ELEV  *'  90  degrees). 

PST  is  the  time-dependent  reference  phase  shitt  per  scan  (in  degrees) 
for  a  given  incoming  signal  condition. 

AN01S  is  the  relative  RMS  value  of  the  noise  signal. 

RN  is  the  starting  number  for  the  pseudo- random  number  generator. 

(Set  RN-0  for  normal  operation.) 

SECT  Is  the  sector  center  (0  £  SECT  <■  360  degrees)  for  a  given  data  run. 
SCANT  is  the  time  taken  to  make  one  complete  scan  (in  seconds)  for  a 
given  data  run. 

I SCAN  is  the  number  of  scans  in  a  given  data  run. 

ISAM  is  the  number  of  sample  points  in  each  scan  for  a  given  data  run. 
SAMI  is  the  azimuthal  angle  between  the  boresight  setting  for  two  conse¬ 
cutive  sample  points  in  a  scan  (in  degrees,  0  SAMI  *  (ISCAN  -1)<360 
degrees)  for  a  given  data  run. 

For  the  input  formats  of  these  input  variables  see  the  program  forward 
of  the  f Emulation  program  GCAAS. 
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The  Basic  Operation  ot  the  Program- 

Reduced  Program  Listing  Number  One 

The  first  necessary  operation,  after  defining  the  natural  constants 

(VEC,  PI,  TWOP1  and  CONRAD)  and  setting  the  data  run  reference  number 

IRUN  to  zero,  is  to  input  the  program  operation  variables  (NRUN,  1TAPE 

and  MODE)  Once  ITAPE  and  MODE  are  known,  then  the  output  data  tape 

can  be  properlv  positioned  and  SETNUM  can  be  determined  for  the  present 

execution  of  the  program. 

After  the  output  tape  has  been  positioned  and  SETNUM  has  been  de¬ 
termined,  the  calculations  for  the  first  data  run  are  begun  by  first 
incrementing  the  data  run  reference  number  by  one,  1RUN“IRUN  +  1.  It 
should  be  noted  at  the  outset  that  the  antenna  array  information  (NOA, 

A,  D,  NAS,  CPHI) ,  the  general  signal  information  (NWAVE,  FREQC ,  FREQM, 
PCAM) ,  the  noise  signal  information  (ANOIS,  RN)  and  the  scanning  opera¬ 
tion  information  (SECT,  SCANT,  ISCAN,  ISAM,  SAMI)  are  all  input  once  for 
the  data  run, and  the  incoming  signal  conditions  (WAMP,  WPHI,  AZM,  ELEV, 
PST)  are  input  NWAVE  times  for  the  data  run.  It  is  necessary  to  re¬ 
supply  all  the  input  information  for  each  data  run  in  the  data  set,  even 
if  certain  input  variables  are  to  be  the  same  for  different  data  runs 
in  the  data  sec. 

The  sequence  of  operations  for  the  data  set,  after  incrementing 
IRUN  by  one  is : 

1)  Input  the  antenna  and  general  signal  information. 

2)  Calculate  the  phase  delays  for  the  delay  lines  . 

3)  Input  one  incoming  signal  condition. 

4)  Calculate  the  antenna  element  voltages  for  the  given  signal  conditions- 

5)  Return  to  3  until  the  calculations  for  all  of  the  incoming  signal 


21 


conditions  have  been  completed 

6)  Input  the  noise  signal  and  scanning  operation  information. 

7)  Define  the  scan  number  IS  and  initialize  the  simulated  time  for  the 
scan. 

8)  Define  the  sample  point  number  IP  and  determine  the  boresight  setting 
for  the  sample  point. 

9)  Increment  the  simulated  time  by  the  time  difference  between  sample 
points  and  determine  the  amplitude  modulation  factor  from  the  sim¬ 
ulated  time. 

10)  Perform  the  time- dependent  phase  shift  calculations  for  the  sample 
point. 

11)  Calculate  the  scanning  bank  voltages  for  the  sample  points. 

12)  Calculate  the  output  data  for  the  sample  point. 

13)  Return  to  8  until  all  of  the  calculations  for  all  of  the  sample 
points  in  the  scan  have  been  completed. 

14)  Return  to  7  until  all  of  the  calculations  for  all  of  the  scans  in 
the  data  run  have  been  completed. 

If  IRUN  f  NRUN,  after  the  calculations  for  the  data  run  have  been 
completed  increment  IRUN  by  one  and  begin  the  calculation  for  the  next 
data  run.  If  IRUN  »  NRUN,  the  execution  of  the  program  is  terminated. 

During  the  operation  of  the  program,  for  a  data  run,  several 
expressions  appear  several  times  or  in  operations  that  are  repeated  many 
times;  in  an  attempt  to  reduce  the  execution  time  for  the  program  these 
expressions  have  been  reduced  to  constants  for  the  data  run  and  incoming 
signal  conditions-  The  data  run  constants  are  defined  as: 
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N0A4  -  NOA/4 
NAS 2  =  NAS/2 

ANOA  -  360/N0A  The  spacing  between  array  elements 

ANOAP  =  ANOA/2  One  hall  the  spacing  between  array  elements 

WL  *=  VEC/FREQC  The  wavelength  of  the  incoming  signals 

CON  ST  1  =■  TWO  PI  *  D/WL 

C0NST2  =  TWOP1  *  A/WL 

CONW  ■  TWOPI  *  FREQW  *  1000  The  angular  frequency  of  amplitude 

modulation 

ANOISJ  =  ANOIS  *  .707  The  relative  RMS  amplitude  of  the 

real  and  imaginary  noise  signals 

WIDTH  -  (ISAM-1)  *  SAMI  The  width  of  the  sampling  window 

s:  r  -  SECT  -  WIDTH/2  The  first  boresight  setting  in  the 

scan 

PELT  -  SCANT  *  SAMI/360  The  time  interval  between  two  con¬ 

secutive  sample  points . 

The  incoming  signal  constants,  used  in  the  antenna  element  voltage 
calculations,  are  defined  as: 

WPH  -  WPHI/ CONRAD  The  reference  phase  in  radians 

COSELV  -  COS (ELEV/CONRAD) 

CQSAMP  -  COSELV  *  WAMP 
C0NE1  ■  C0NST1  *  COSELV 
CONE 2  -  C0NST2  *  COSELV  . 

Both  the  boresight  settings  and  the  simulated  time  are  calculated 
for  the  sample  point  from  the  scan  and  sample  point  reference  number. 

The  boresight  setting  can  be  found  from  IP  alone,  since  the  boresight 
settings  will  always  be  the  same  for  each  scan  in  a  data  run.  The  bore¬ 
sight  setting  for  sample  point  IP  is  determined  as: 

CENT  -  (IP-1)  *  SAMI  +  SIDE. 
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It  should  be  noted  that  It  may  be  necessary  to  correct  CENT  to  lie 
between  the  limits  0  -  CENT  360  degrees. 

The  simulated  time  is  computed  by  Initializing  the  simulated  time 
at  the  beginning  of  the  scan. 

T  =  (IS  -  1)  *  SCANT  -  DELT,  and  then  increment  the  simulated 

time  by  the  time  interval  between  sample  points,  T  =  T  +  DELT,  for  each 

sample  point.  It  is  interesting  to  note  that  this  method  of  computation 

of  simulated  time  is  equivalent  to: 

T  *  (IS  -  1)  *  SCANT  +  (IP  -  1)  *  DELT. 

Once  the  simulated  time  is  known  for  the  sample  point,  the  amplitude 

modulation  factor  is  determined  as: 

AMS  -  (1  -  PCAM)  *  PCAM  *  COS  (CONW  *  T) 

which  is  Equation  (14),"^ 

Phase  Delay  Calculations- 

Reduced  Program  Listing  Number  Two 

The  first  step  in  the  phase  delay  calculation  is  to  reduce  redundant 
calculation  by  reducing  the  two  terms  in  Equation  (17) ,  which  are  the 
same  for  all  the  scanning  elements,  to  two  constants.  These  two  con¬ 
stants  are  defined  as: 

CONCHI  -  COS  (CPHI /CONRAD)  *  C0NST2 

COSM  =  COS  (((NAS 2  -  1)  *  ANOA  +  AN0A2) /CONRAD) . 

Equation  (17)  is  then  implemented,  for  scanning  element  ISE,  as: 

PDELAY  -  CONCHI  *  (C0S(((ISE  -  1)  *  ANOA  +  AN0A2) /CONRAD)  -  COSM). 

The  phase  delay  information  is  then  stored  as: 

COSP  (ISE)  =  COS (PDELAY) 

SINP  (ISE)  -  SIN (PDELAY) , 


I  Equation  numbers  refer  to  equations  in  Chapter  3. 
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since  the  sine  and  the  cosine  of  the  phase  delays  will  be  the  only 
quantities  used  in  further  calculations. 

Antenna  Element  Voltage  Calculatlon9- 

Reduced  Program  Listing  Number  Three 

After  the  conditions  have  been  read  and  the  repetitive  termB  used 
in  the  antenna  element  voltage  calculation  have  been  reduced  to  con¬ 
stants  for  the  incoming  signal  IWAVE,  the  antenna  element  voltages  ex¬ 
cited  by  the  given  signal  conditions  are  determined  in  the  following 
manner : 

1)  Determine  the  element  K  to  the  left  of  the  azimuthal  angle  of 
arrival  as  K  =  (AZM/ANOA)  +  1. 

2)  Determine  the  spacing  between  the  azimuthal  angle  of  arrival 

and  the  first  element  to  the  left  as  DELK  ■  AZM-ANOA  *  (K-l) 

and  the  first  element  to  the  right  as  DELKP  *  ANOA-DELK. 

3)  For  the  elements  to  the  left  of  the  azimuth  IDIR  ■  1. 

4)  Set  IAW  =  1  for  the  first  element  to  the  left  of  the  signal 
azimuth. 

5)  Find  the  proper  antenna  reference  number  as  I ANT  »  K  +  1-IAW. 

6)  Determine  the  spacing  between  the  azimuth  and  the  element  as 
THE  =  ANOA  *  (l-IAW)-DELK. 

7)  Implement  Equations  (1)  and  (2)  to  find  the  voltage  excited 
by  the  signal. 

8)  Implement  Equations  (3)  and  (4)  to  find  the  total  element 
complex  voltage. 

9)  For  the  next  element  to  the  left  IAW  ■  IAW  +  1. 

10)  Return  to  5  if  IAW  s  N0A4. 
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Note  that  voltages  will  only  be  excited  in  antenna  elements  up 
to  one  iourth  ot  the  distance  around  the  array  from  the 
azimuthal  angle  ot  arrival 

11)  For  the  elements  to  the  right  of  the  azimuth  ID1R  •  2. 

12)  Set  1AW  -  1  ror  tne  first  element  to  the  right  of  the  signal 
azimuth 

13)  Find  the  proper  antenna  reference  number  as  1ANT  •  K  t  1AW. 

14)  Determine  the  spacing  between  the  azimuth  and  the  element  as 
THE  -  ANOA  *  UAW  —  1)  r  DELKP. 

15)  Implement  Equations  (.1)  and  f2)  to  find  the  voltage  excited  by 
the  signal. 

16)  Implement  Equations  C3)  and  (.4)  to  find  the  total  element 
complex  voltage. 

17)  For  the  next  element  to  the  right  IAW  ■  IAW  +■  1. 

18)  Return  to  13  it  IAW  ^  NGa4. 

It  should  be  noted  that  It  may  be  necessary  to  correct  IANT  and  THE 
to  fit  the  limits  i  IaNT  i  NOA  and  j THE  |  ^90  degrees,  before  being 
used  in  any  calculations. 

Equation  CD  and  (2)  are  implemented  in  the  program  as: 

AMP  -  COSAMP  *  SIN  (C0NE1  *  COSTHE) 
and 

PHI  =  CONE 2  *  COSTHE  -r  WPH 

where  COSTHE  is  defined  for  the  given  element  IaNT  as: 

COSTHE  -  COS (IHE/ CONRAD) 

and  where  COSAMP,  C0NE1,  C0NE2  and  WPH  are  defined,  in  reduced  listing 
number  one,  for  the  given  incoming  signal,  IWaVE. 
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After  Equations  (1)  and  (2)  have  been  implemented,  Equations  (3) 
and  (4)  can  then  be  implemented,  for  antenna  element  1ANT  and  incoming 
signal  IWAVE ,  as: 

ARE  (1ANT)  -  ARE  (IANT)  +  AMP  *  COS'(PHI) 

and 

AIM(IANT)  -  AIM  (IANT)  t  AMP  *  SIN  (PHI). 

ARE(IANT)  and  AIM ( IANT ) ,  the  antenna  element  voltage  variables,  were 
zeroed  for  all  NOA  elements  before  the  calculations  for  the  first  in¬ 
coming  signal  conditions  were  begun. 

For  the  implementation  of  the  time-dependent  phase  shift,  it  will 
be  necessary  to  store  certain  intermediate  values  from  the  antenna 
element  voltage  calculation.  These  values  are  stored  as: 

WPHH( IWAVE)  =  0.0,  the  old  reference  phase  shift  for  Bignal  IWAVE; 

PSTEP (IWAVE)  =  PST/CONRAD,  the  reference  phase  step  per  scan,  for  signal 
IWAVE; 

PHIH( IANT, IWAVE)  -  PHI,  the  base  phase,  for  element  IANT  and  signal  IWAVE; 
AMPH(IANT,  IWAVE)  =  AMP,  the  base  amplitude,  for  element  IANT  and 
signal  IWAVE,  which  is  also  zeroed  before  calculations  for  the  first  in¬ 
coming  signal  conditions  were  begun. 

Time  Dependent  Phase  Shift  Calculations- 

Reduced  Program  Listing  Number  Four 

The  time-dependent  phase  shift  calculations  are  carried  out  for 

each  sample  point  for  those  incoming  signals  which  have  a  time  dependent 

reference  phase  shift,  i.e.  PSTEP(IWAVE)  jrf  0.  It  will  also  only  be 

necessary  to  carry  out  the  calculations  for  those  elements  which  have 

a  non-zero  base  amplitude  for  the  given  signal,  i.e.  AMPH (IANT, IWAVE) 


The  first  step  in  the  calculations,  for  a  given  signal  1WAVL ,  is 
to  determine  the  reference  phase  shift,  for  the  present  sample  point  at 


Simulated  time  T,  as:  WPPH  =  PSTEP(IWAVE)  *  T/SCANT.  The  reference  phases 
for  Che  past  and  present  sample  points,  Equation-  (.7)  and  (8),  are  then 
obtained  by  adding  the  element  base  reference  phase: 

PAl  -  PH1H  UANT,  IWAVE)  r  WPHH(IWAVE) 

and 

PA 2  =  PH1H  (IANT,  1WAVE)  -r  WPHH  . 

Equations  (5)  and  (6)  can  then  be  Implemented  aB : 

ARE  ( IANT )  =  ARE(IANT)  i-  AMPH(IANT,  IWAVE)  *  (COS (PA2)-C0S (PAl) ) 

and 

AIM  (IANT)  =  AIM(IANT)  i-  AMPH(1ANT,  IWAVE)  *  (SIN (PA2)-S1N (PAl) ) . 

The  last  step  in  the  calculation,  for  the  given  signal,  is  to  store 
the  present  reference  phase  shift  in  the  variable  used  as  the  old 
reference  phase  shift,  i.e,  WPHH(IWAVE)  =  WPPH. 

It  should  be  noted  that  these  calculations  are  bypassed  for  T  ■  0, 
because  the  phase  shift  is  zero. 

Scanning  Element  and  Scanning  Bank  Voltage  Calculations- 

Reduced  Program  Listing  Number  Five 

Once  CENT,  T,  and  AMS  have  been  determined  for  the  sample  point 
and  the  time  dependent  phase  shift  calculations  have  been  completed, 
the  scanning  bank  voltages  are  calculated  in  the  following  manner : 

1)  Determine  the  antenna  element  KC  to  the  left  of  the  first  ele¬ 
ment  in  the  left  scanning  bank  as  KC  ■  (CENT-AN0A2) /’ANOA  +  1. 

2)  Determine  the  spacing  between  a  scanning  element  and  the  first 
antenna  element  to  the  left  as  CDEL  -  (CENT-AN0A2)- (ANOA  *  (KC-1)) 
and  the  first  antenna  element  to  the  right  as  CDELP»ANOA-CDEL . 


L 
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3)  Determine  the  linear  displacement  factor  for  antenna  elements 
to  the  lef.  as  CDF.LAP  -  CDELP/ANOA  and  to  the  right  as  CDELA 


4) 

5) 

6) 
7) 


8) 

9) 

10) 

ID 

12) 

13) 


14) 

15) 

16) 
17) 


CDEL/ANOA. 

For  the  scanning  elements  in  the  left  bank,  IDIR  ■  1. 

Zero  the  bank  voltage  accumulator  SUMI  ■  SUMR  ■  0. 

Set  ISE  =  1  for  the  first  element  in  the  left  bank. 

Find  the  proper  antenna  element  reference  number  for  the  ele¬ 
ment  to  the  left  of  ISE  as  IL  ■  (KC-ISE  +  1)  and  the  element 
to  the  right  as  IR  =  (IL  +  1) . 

Determine  the  noise  voltage  for  tne  element  to  the  left  as 
ANREL,  ANIML,  and  the  element  to  the  right  as  ANRER,  ANIMR. 
Implement  Equations  (15)  and  (16)  to  find  the  scanning  element 
voltage. 

Implement  Equations  (18)  and  (19)  to  find  the  scanning  bank 
voltage. 

For  the  next  scanning  element  in  the  left  bank  ISE  ■  (ISE  +  1). 
Return  to  7  if  ISE  <;  NAS2. 

Store  the  left  scanning  bank  voltage  as: 

SUMRE (IDIR)  -  SUMR 
SUMIM(IDIR)  -  SUMI 

For  the  scanning  elements  in  the  right  bank  IDIR  ■  2. 

Zero  the  bank  voltage  accumulator  SUMI  ■  SUMR  ■  0. 

Set  ISE  ■  1  for  the  first  element  in  the  right  bank. 

Find  the  proper  antenna  element  reference  number  for  the  ele¬ 
ment  to  the  left  of  ISE  as  IL  =  (KC  +  ISE)  and  the  element  to 
the  right  as  IR  ■  (IL  +  1). 
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18)  Determine  the  noise  voltage  tor  the  element  to  the  left  as 
ANREL,  AN1ML  and  the  element  to  the  right  as  ANRER,  AN  1 MR . 

19)  Implement  Equations  (15)  and  (16)  to  find  the  scanning  element 
voltage 

20)  Implement  Equations  (.18)  and  (19)  -o  i  ind  the  scanning  bank 
voltage . 

21)  For  the  next  scanning  element  in  the  right  bank  ISE  -  (1SE  +  1). 

22)  Return  to  17  if  ISE  r,  NaS2. 

23)  Store  the  right  scanning  bank  volLage  as: 

SUMRE  (1D1R)  *  SUMR 

SUM1M  (1D1R)  -  SUM! 

Again  it  should  be  noted  that  It  may  be  necessary  to  correct  1R 
and  IL  to  fall  in  the  limits  1  co  NOA  before  they  are  used  in  any  cal¬ 
culations  . 

Equations  (15)  and  (16)  are  implemented  as: 

ERE  -  (AMS  *  ARE (IL)  +  ANREL)  *  CDELAP  +  (AMS  *  ARE ( IR)  +  ANRER)  *  CDELA 

SIM  -  (AMS  *  AIM(IL)  +  ANRML)  *  CDELAP  +  .(AMS  *  AIM(IR)  +  ANIMR)  *  CDELA 

Equations  (18)  and  (19)  are  implemented  as: 

SUMI  -  SUM1  -t-  SIM  *  COSP(ISE)  -  SRE  *  SINP(ISE) 

SUMR  -  SUMR  +  SIM  *  SINP(ISE)  +  SRE  *  COSP(ISE) 

where  SINP(ISE)  and  COSP(ISE)  are  the  stored  sine  and  cosine  of  the 

phase  delay  for  element  ISE. 

Determination  of  the  Noise  Signal  Voltage- 

Reduced  Program  Listing  Number  Six  and  Seven 

The  noise  voltages  are  determined  by  making  two  successive  calls  to 

the  computer  library  pseudo-random  number  generator  (CALL  RND(RN,  RNDP)) 

and  multiplying  the  supplied  real  number  RNDP  by  the  relative  noise 
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voltage  factor  AN01SJ,  where  RNDP  Is  from  a  GauBSion  distribution  and 
RN’DP  i  5.0.  It  is  necessary  to  determine  the  noise  voltages  for  the 
antenna  elements  to  the  right  and  left  of  the  scanning  elements  in 
different  manners  for  the  left  and  right  bank. 

For  the  left  bank,  since  the  order  in  which  calculations  are  made 
for  the  scanning  elements  proceeds  from  right  to  left,  the  noise  voltages 
ANREL,  ANIML  for  element  (ISE-1)  become  the  noise  voltages  ANRER,  ANIMR; 
for  element  ISE  and  ANREL,  ANIML  is  determined  by  consecutive  calls  to 
the  random  number  generator.  For  (ISE*  1),  it  is  necessary  to  make 
four  calls  to  the  random  number  generator  to  determine  ANRER,  ANIMR 
and  ANREL,  ANIML;  it  will  also  be  necessary  to  store  ANRER,  ANIMR  as 
ANRERF,  ANIMRF  since  the  antenna  element  to  the  right  of  the  first  ele¬ 
ment  in  the  left  bank  is  also  the  antenna  element  to  the  left  of  the 
first  element  in  the  riglit  bank. 

For  the  right  scanning  bank,  since  the  order  in  which  calculations 
are  made  for  the  scanning  elements  proceeds  from  the  left  to  the  right, 
the  noise  voltages  ANRER,  ANIMR  for  element  (ISE-1)  become  the  noise 
voltages  ANREL.,  ANIML  for  element  ISE  and  ANRER,  ANIMR  is  determined  by 
consecutive  calls  to  the  random  number  generator.  For  ISE  ■  1,  the  noise 
voltage  ANREL,  ANIML  is  set  equal  to  ANRERF,  ANIMRF  and  the  noise  voltage 
ANRER,  ANIMR  is  determined  by  the  random  number  generator. 

Calculation  of  System  Outputs- 

Reduced  Program  Listing  Number  Eight 

After  the  bank  voltages  have  been  determined,  the  output  data  are 
determined  for  the  sample  point. 

The  sum  data  are  found  as  the  magnitude  of  the  complex  sum  of  the 
right  and  left  bank  voltages,  by  implementing  Equation  (20)  as: 
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SUMR  -  SUMRE  (2)  i-  SUMRE(l) 

SUMI  -  SUMIM(2)  +  SUMIM(l) 

DATAST  «>  SQRT (SUMR  **  2  +  SUMI  **  2). 

The  difference  damn  la  found  as  the  magnitude  of  the  complex 
difference  between  the  right  and  left  bank  voltage,  by  implementing 
Equation  (21)  aa : 

SUMR  -  SUMRE (2)  -  SUMRE(l) 

SUMI  =  SUMIM(2)  -  SUMIM(l) 

DATADT  =  SQRT (SUMR  **  2  +  SUMI  **  2). 

The  differential  phase  dattan is  found  as  the  difference  between  the 
phase  of  the  right  bank  voltage  and  the  phase  nf  the  left  bank  voltage, 
by  implementing  Equation  (22)  as 

ANGR  -  ATAN2  (SUMIM(2) ,  SUMRE (1)) 

ANGL  -  ATAN2  (SUM1M(1),  SUMRE(2)) 

DATP  -  ANGR  -  ANGL  . 

where  ATAN2  (IM,  RE)  is  defined  over  the  range  -180  '  ATAN2  (IM,  RE) 

<  +180.  It  may  be  necessary  to  correct  DATP  to  fall  between  the  limits 


-180  <  DATP  *  +180. 
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V.  OUTPUT  DATA 

The  simulation  program  produces  both  a  hard  copy  and  magnetic  tape 
output . 

On  the  first  page  of  the  hard  copy,  for  each  data  run,  the  input 
conditions  for  the  simulation  are  defined  along  with  the  data  set  and 
data  run  numbers  which  are  used  in  locating  the  data  run  on  the  magnetic 
tape.  Succeeding  pages  contain  the  output  data  for  each  scan  in  the  data 
run. 

After  the  calculations  for  each  sample  point  have  been  completed, 
as  discussed  in  the  last  chapter,  the  output  data  for  sample  point  IP  are 
stored  as: 

DATAC  (IP)  =  CENT 
DATAS  (IP)  -  DAtEAST 
DATAD  (IP)  -  DATADT 
DATAP  (IP)  =  DATP  . 

As  the  sample  point  data  are  being  stored,  the  sum  and  difference  data 
are  tested  to  determine  their  maximum  values  for  the  scan;  these  maximum 
values  are  then  stored  as  DATASM  and  DATADM, respectively .  After  all  of 
the  calculations  for  the  scan  have  been  completed,  the  output  data  for 
the  scan  are  then  printed  on  hard  copy. 

Along  with  the  numerical  values  of  the  output  data,  a  fifty-one- 
point  incremental  line  printer  plot  is  also  displayed  on  the  hard  copy. 
The  sum  and  difference  data  are  normalized  by  DATASM  and  DATADM  respec¬ 
tively,  and  then  plotted  with  the  maximum  values  as  the  right-hand  side 
of  the  plot.  The  differential  phase  data  are  plotted  between  the  limits 
-180  degrees  on  the  left  and  +180  degrees  on  the  right  hand  side  of  the 
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plot.  The  symbols  used  in  the  plot  to  represent  the  output  data  are: 

*  Sum 

+  Difference 

X  Differential  Phase. 

Appendix  C  contains  the  hard  copy  output  for  a  sample  execution  of  the 
program;  the  necessary  input  cards  for  the  execution  are  also  shown. 

The  organization  of  the  output  data  tape  is  such  that  each  data 
set  will  comprise  one  tape  file.  After  the  file  for  the  last  data  set 
on  the  tape, there  is  a  terminator  file  which  is  used  in  the  tape  han¬ 
dling  procedure.  The  data  set  file  consists  of  a  series  of  logical 
records,  one  H  (or  input)  record  for  each  data  run  followed  by  one  S 
(or  output)  record  for  each  scan  in  the  data  run  and  a  terminator 
file  which  consists  of  one  E  (or  end)  record.  The  structure  of  the  H, 

S  and  E  records  is  given  in  Tables  1,  2,  and  3,  respectively. 

The  simulation  program,  as  was  noted  in  the  last  chapter,  has  two 
modes  for  tape  handling  operations;  the  first  mode  is  used  for  new  tapes 
and  the  second  mode  is  used  for  old  tapes  to  which  an  additional  file 
is  to  be  added.  In  the  first  mode  it  is  only  necessary  to  open  the 
tape  and  position  it  to  its  beginning  and  then  the  first  data  set  file 
is  written  on  the  tape.  The  second  mode  is  used  to  position  the  tape 
to  the  end  of  the  last  data  set  file  on  the  tape;  this  is  done  by 
searching  the  tape  for  the  terminator  file  and  then  writing  the  new 
data  set  file  over  it.  In  both  modes  one  and  two,  after  the  new  data 
set  file  has  been  added  to  the  tape  the  terminator  file  is  then  added 
to  the  tape.  The  only  information  in  the  terminator  file  is  the  set 
number  of  the  last  data  set  file  which  is  used  in  mode  two  to  determine 
the  set  number  for  the  next  data  set  file  to  be  added  to  the  tape. 
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Table  1 

The  input  Information  for  the  data  run  is  organized  In  the  H  record 


In  the  following  manner: 


AREA(  1) 

4H 

AREA(  2) 

SETNUM 

AREA(  3) 

IRUN 

AREA (  4) 

NRUN 

area(  5) 

NOA 

AREA(  6) 

A 

AREA(  7) 

S3 

D 

AREA(  8) 

NAS 

AREA(  9) 

CPHI 

AREA (10) 

FREQC 

AREA  (ID 

B 

FREQM 

AREA (12) 

= 

PCAM 

AREA (13) 

NWAVE 

AREA (14) 

ANOIS 

AREA (13) 

RN 

AREA (16) 

SECT 

AREA (17) 

ISCAN 

AREA (18) 

SCANT 

AREA (19) 

ISAM 

AREA (20) 

SAMI 

AREA (21) 

WAMP 

AREA(22) 

WPHI 

AREA(23) 

AZM 

AREA(24) 

ELEV 

AREA(25) 

PST 

AREA(26) 

WAMP 

AREA (27) 

WPHI 

AREA (28) 

AZM 

AREA (29) 

ELEV 

AREA (30) 

PST 

AREA (31) 

S 

AREA(  ) 

AREA(  ) 

AREA(  ) 

-  Set  number 

-  Run  number 

-  Number  of  runs  in  the  set 

-  Number  of  antenna  elements  in  the  array 

-  Radius  of  the  reflecting  cylinder 

-  Cylinder  to  element  distance 

-  Number  of  scanning  elements 

-  Cophasel  angle 

-  Carrier  frequency 

-  Modulation  frequency 

-  Modulation  factor 

-  Number  of  incoming  signal  conditions 

-  Relative  noiBe  level 

-  Starting  Number  for  pseudo  randum  number  gen. 

-  Sector  center 

-  Number  of  scans  in  the  run 

-  Time  to  complete  on  scan 

-  Number  of  sample  points  in  each  scan 

-  Spacing  between  sample  points _ 

-  Relative  signal  amplitude 

-  Reference  phase  For 

-  Azimuthal  angle  of  arrival  IWAVE  -  1 

-  Elevation  angle  of  arrival 

-  Time  Dept,  phase  shift _ 


For 

IWAVE  -  2 


To  access  the  information  for  signal  conditions  IWAVE: 


WAMP 

WPHI 

AZM 

ELEV 

PSR 


AREA (21  +  (IWAVE-1)  *  5) 
AREA(22  +  (IWAVE-1)  *  5) 
AREA(23  +  (IWAVE-1)  *  5) 
AREA (24  +  (IWAVE-1)  *  5) 
AREA (25  +  (IWAVE-1)  *  5) 
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Table  II 


The  outpuc  information  for  a  scan  is  organized  in  the  S  record  in 
the  following  manner: 


AREA(  1) 
AREA(  2) 
AREA(  3) 
AREAC  4) 
AREA(  5) 
AREA(  6) 
AREA(  7) 
AREAC  8) 
AREAC  9) 
AREA (10) 
AREA(ll) 
AREA(12) 
AREA (13) 
AREA (14) 
AREA (15) 
AREA (16) 
AREA (17) 
AREA(18) 
AREA(19) 
AREA (20) 
AREA (21) 
AREA(22) 
AREA(23) 
AREA (24) 
AREA(  ) 
AREA(  ) 
AREA 


4H  S 

SETNUM  -  Set  number 

IRUN  -  Run  number 

NRUN  -  Number  of  run9  in  the  scan 

IS  --  Scan  number 

DATASM  -  Maximum  value  of  the  sum  data 
DATADM  -  Maximum  value  of  the  difference  data 


Blank 


CENT 

DATAST 

DATADT 

DATP 

CENT 

DATAST 

DATADT 

DATP 

CENT 

DATAST 

DATADT 

DATP 


Boresight  setting 

Sum  data 

For 

Difference  data 
Differential  phase  data 

IP-1 

Boresight  setting 

Sum  data 

For 

Difference  data 
Differential  phase  data 

IP-2 

For 

IP-3 

To  access  the  output  data  for  sample  point  IP: 

CENT  -  AREA (11  +  (IP-1)  *  4) 

DATAST  -  AREA (12  +  (IP-1)  *  4) 

DATADT  -  AREA(13  +  (IP-1)  *  4) 

DATAP  -  AREA (1 4  +  (IP-1)  *  4) 


Table  III 


The  E  record  is  organized  in  the  following  manner: 
AREA(  1)  =  4H  E 

AREA(  2)  ■  SETNUM  -  Set  number  of  last  data  set  file 
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The  method  chat  should  be  used  to  locate  a  given  data  run,  defined 
by  SETNUM,  IRUN,  is  to  test  the  first  record  in  each  file  until 
AREA(2)  -  SETNUM  and  then  check  each  record  in  the  file  until  AREA(l)  ■ 
4H _ H  and  AREA(3)  -  IRUN.  Appendix  D  contains  a  segment  of  a  pro¬ 

gram  which  can  be  used  to  properly  locate  a  given  data  run. 

Appendix  E  contains  the  listing  of  program  STRO,  standing  for 
Simulation  Tape  Read  Out,  which  is  used  to  obtain  a  hard  copy  print 
out  of  all  the  information  on  the  output  tape. 
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c 

(-10.? 

i  m  c  a  n  r 

-  1 

.up  [  pp  f  *,l  v,\N  1  N  SFC'lNOC 

R  A  III  R 


o  o  o  o  o 
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c 

c 

c 

c 

C 

C 

C 

C 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

t: 

c 

c 

c 

c 


c 


c 

c 

c 


OTP 

Oil 

0 

c 

c 


c 

c 

r 


(  I  S  )  I  SCO 

(IS)  ISA* 

(  '  1  n  .  ?  )  s  a  m  ! 


NUMBER  m-  SI  ANS  IN  THF  11ATA  RL  N 
NiimRFP  Oi-  S'MPI.E  POINTS  IM  EACH  SCAN 
ANRLF  BtHc.N  SAMP!  F  POINTS  IT  Q-GREES 


P  T  0  G  w  A  N’  I  IT’ ITS 
I  MPtiT : 

MAXIPUu  NUMBER  I'.f  SIGNAL  ('PM  I  T  I  ()NS  S 

A  7  Ml  I  T  H  A  I  A.'GlrS  SHOULD  RE  EXPRESSED  RPTWFPN  THF  L  1  M  I  T  S  OT 
O.n  ofgr-FS  AM*  350.99-  DPGRffs 

0  J  T P l;  T  : 


MAXIMUM 

\!  U  M  P  R  R 

mF 

SAMP|  E  FC  I  ' Ts  /  SCAN 

inn 

0  JFR  a  T  . 
M A  x I T  UM 

on  : 

number 

:F 

ANTENNA  El  1  MFNTS 

36n 

MAX  I  MR- 

JUMPER 

F 

scanning  a- tfmna 

ion  (mlst  he  an  bvt-T.  t c  .  > 

#**##*■#  u  t' #«#*#####«  ###*######■»####*****##•  «»■>** 

T4PC  M  M'M.  I  T'C- 

THF  0UTP"T  DATA  U  1 1  L  BF  STORF'l  0M  THF  TAPP  M0UTTE'.  ON  l)«]Vt  |T»r--t 

THF  r |  p  ■  aTA  TAPE,  IT  ANY,  SfOUL'I  HP  MOUNTED  ON  DRIVE  lUh- 
T  h  f  U-’F  Will  Hi-  AiVANCED  TO  ThR  FNn  OT  THf  CLC  DATA  A.'.o  T  h6  t  F « 
Data  Wrll  FF  ST  f.RFI'  BEGINNING  AT  that  POINT  ON  THF  TApF 


***#«*»»******»*«» 


i*«*«**#**i 


o  IMF' SI  0“-  A  PRAY' 


n  i  mf' s  i  n“- 

I )  1 m  F ' S i ON 
D  i  mf' s  i  ri'i 
0  I  MF'  S  O'l 
D  I  mp  s  •  ns 
n  I  m  c  in  s  n ' 
0  I  M F '  S  ON 
D I  mc'.s  i  n 


mfP(SO)  ,SINP<5U) 

ARF  (3*0  )  .  A  I  M(  3 (S 0  ) 
SIIMRF  I  2  )  ,  SOM  T  M  <  2  ) 

r  atas  ( i  nn ) ,  tiat at  ( i  n o 

WPHH(t-)  ,  PSTFp  (  5  ) 

PH  I  H (  sen ,  s  )  ,  A  M  PF <  V" 
A |  I.NF(51',AMARK(11' 

A  R F  A  (  a  1  D  )  ,  DF  A  N  (  6  5  0 


i  .  i 


A  T  A  P  (  1  o  (1  >  ,  D  A  T  A  0  ( 1  0  0  ) 


5  ) 


H  H  =  4  F  h 
S  h  =  4 1 

F  h  =  4 ^ 

I  \i  I  T  1  A  /  F  II  N  F  F  P  Ti  T  F  R 

00  or  n  i  P:i  ,  S] 

A  .  J  M-  (  ;  *-  1  :  P 

DO  9l>  1  I  °r1  .1  1 

A  M A R -  I  1  >  ) -  27 


ii-p  r  p 


.  T  1 10  A  l. 


CNF  r ANTS 


VrCsOOO . h 
31 =7 ,1 41 S  9?  *  p 
T  a, no  =  I  »P 
COMO  n-TAfl/’,  )F 


r o.  -|.,n  v  i;p,  rai|(t  Ini-mo  ai!  n 


ooo  oorj  o  o  o  ei 


HUMrO 

R  -  a  n  go  ofl  . r 'Wl'\ ,  1  T  a^E  ,  mode 
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C 

C  0jem  T  ■  p*-  0\  I)R  VE  I  T  a  PE 

C 

D  ■"  A  S  =  6  n  o 

CALL  SETHNI  I T A Pt  ,  1  ) 

C  4  L  (  fvRLKI  I  T  A  Pf-  ,  1  ) 

CALI  RL  KLN(  ITAPE.0.0) 

C 

G1  Tn  ( Q?4, 914  )  , MODE 
C 

C  TAPP  M'  0t  NO  .  ? 

c  a  i v a A'C l-  tape  to  em  or  on:  data 

c 

914  CONTINUE 

915  A  3  S I C  M  6  8  0  0  TO  IASI 
91ft  CALI  O-FN ( PPAN , '  FAS, ]  TAPE  .  0  .  4EIS  !  "'ll  > 

CALL  R  >LR(0FAN. AREA, LAST) 

CALI  Cl  OSE ( OF  AM ) 

Ir  liRMIl  >-FH)  517,916.917 
917  C  *  L  l.  S^TFD(0F  AN) 

GO  TO  91ft 

9 1 P  CAl.i  SFTPK(OFAN) 

C  A I  L  rWEN<nFAN,r FAS, ITAPE.0.4HSIMU) 
SETni  ^=apFa (?) 

GO  Tf’  9^0 

TAPP  NO.  1 

924  C3NTIMMF 

925  CALI  0 - = \ C  i  FAN, -FAS, I  TAPE . 0 . 4MS t ) 

s  =  t  m  i  ■  m  =  n 

930  S?TMI:M  =  SPTNUM*inOO 

8  r  G  T  \  N  •  M f'  OF  0  A  r  A  RUN 
950  I ?UN= I RUN*} 


R-An  A'T-nna  IM-CR^ATION  ANn  GENcRal  SIGNAi  INFORMATION 
C 

R-An  9009, NOA, A . C  » *  AS.CPHI 
W  -  a  n  onin,  Ni  wave  .FRi-OC,FRECM,pcam 
c 

C  DETERMINE  DATA  ‘IN  OPERATION  CONSTANTS 

C 

M  D  A  4  =  N  A  /  4 
A\GA=T^0/NP6 
A  Ml)  A  5s  N  J  A  /? 

NAS?rN/iR/? 

*1 .  S  w r-  r  /  r  k  P  0  c 

C  J f  s  1  1  =  T  .V  n  P  T  •  p  /  L 
C  0 1 ' 9  T  2  =  T  n  0  P  I  *  A  /  w  L 
COMWrTwn-  I  *FRF  Uv  *  1 1.  0  0 
C 

c:  fjoi  r  i  -  T-  PRASE  CPI  AYS  FOR  S  4NM;NKi  AM-NNa  FI.F^F  ,TS 

C 

C  j  Mr  *■  I  -r  lAfpPR  I  /  CGI1  R  A  0  )  «  C  r  N  S  i  2 
CnPMrC  S  (  (  f  N  A  S  2  -  1  )  *ANOA*ANOA5  )  /  '  UNR  A  I)  ) 


o  o  o 
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DO  909  I SF  =  1 , NA- 2 

PDF)  />  Vriy'NTH  I  •  (I  CS(((!SE-l)«<'Nr)A*A\0*?)/CnNRA'))-Cl)SM) 
C  DSP  <  !  'F ) sCCS ( P  6LAY) 

S  l  NP  (  I^F  >*SIN  tPi'El  ay) 

999  CONTINUE 
C 

C  PRINT  r  HECORn  1  KK  IiRM  A  T  I  on 

c 

PRINT  9f)19»SFTNljM.  1  R  U  N  ,  N  W  l  N 

PRINT  9n?P.MOA,  A,  D,  NAS.  CPH  ,  f  RtQC,FREOM  ,  PCAM,N|»AV- 
P  5  I  M  T  9  fl 1 
C 

C  STORE  h  record  information 

c 

AREA (  1 ) =HH 
AREA!  ?  )  =  SE  TNIJM 
A  RE  a  (  .t»  =  IRUN 
AREA)  4):NE0N 
AREA(  SjsNlOA 
A  R  F  A (  6)=A 

A  RE  A (  7)=D 
AREA!  fl  )  =  N  A  $ 

AREA)  9  )  =  C  P  H I 

area ( in )=FPFPO 
ARt A ( 1 1  )  =FREOM 
AREA  < IP  )  =Pf  AM 
A  R  E  A  (  1  .7  )  =  N  k  A  V  E 
K<K=?1 
C 

C  7  ERO  ANTENNA  ELEMENT  ARRAYS 

C 

m  1  non  I  AMT  si , MCA 
AREf IA JTlid.O 
A  f  M  f  |  A  f  1 T  )  s  0 .  n 
DD  iron  IWAVEsl.NWAVE 
A  4  PEI  1  I  iN-i'WAVE'sP.O 

1  n 0 n  cD'mtimiif 
c 

I  <1  A  VE  =  1 
C 

c  R  =  AO  A  D  PRINT  SIGNAL  CONT  I  T  i  ON 

C 

1  nm  R  r  A  n  9ni1fRAMP,^PHItA7M,EI-E'/.PST 

P9  I  MT  onVP,?WAV'-,WAMP,WPHI*A/M,Fi  EV , PST 

C 

wjw:,PM[/r(MRAi) 

CTRC  V-C  S  (  Et.F  V/CON'RAD) 

CDSavPsO'tSEI.  V  *  w  a  M  P 
CDMP"  =  'n  \flTl  *  EL  v 
CDNC9S  "  S  T  ?  •  0  0 E L 
w j  h  m  1 1  >i  /p  )  sp .  n 
P  S  T  c  -•  (  ^  n  \/  p  )  s  P  S  /C>  NR  AD 

sr-'Di  .. crr.cn  nflrmation 

A  RE  A  t  R'R  )  sW  AM- 
A  }  F  A  (  K  •'  <  *  1  IsEPm: 

A  Rp  4 ( K • *  +  ?  )sA7M 
A  3E  A  I  K  ■  K  +  T ) sP!  t  ' 

A  RE  A ( K  '  <*4 >  =ES I 

R  <  4  s  "  x  ♦  fl 
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r 

C 

C 

r 


c 

c 

r 

c 

c 

c 

1002 


C 

C 

C 

1003 

10  04 

1005 

1006 
1007 

C 

C 

C 

100P 

1009 
10  10 

10  11 
C 
C 
C 

101? 


C 

C 

C 


1013 

1014 


1  016 

c 

c 

r. 


c 

c 

c 


O-TPBM  \j  uT(jN  A  i  L  b  M  P  \  T  TO  T  H  E  L  tr  F  T  {A  T  H  E  S  I  li  N  A  L  A  Z  1  M  U  '  H 

AMO  I  H  I'nlER  -ETpEEN  THF  A7MIHIITH  AM  AO.iAI  fM  ANTENNA  1 1  fc  M  b  N  F  b 

6  :  (  A  7M/AN0A  )  *1 

II  ft  K-A/M-ANOA*  (  '  *1  ) 

D  -  L  K  F>  =  A  N I  A  -  U  E  '  6 


i;i  rf;  p  (■  t  for  I C  In  =  i  ANn  ppm  i  ow  imps? 

no  i  n 1 4  i  n  i r  =  i . 

PREFER"  )PF  R  A  T  I  N  (  OR  ANTENa  HEmENTS  OP  Tf  50  DEG .  FRO*  aZIElTp 

rn  i  n  i  .1  i  A  ^  =  1  ,  'V  I A  4 
G  3  T  ■  1  0  0  3  . 1  0  n  a  )  ,  I  n  1  R 

F  !  A: 0  A  T-NvJA  M  BEE'  F  OP  OEERaT  I  ON  TO  IEaT 
JANTrK+l-.AW 

I  •'  (  U  T)  1004,1004,1005 

I  ANT r MO  A* I  ANT 

T  HE  =  A  NOA  »  < 1  -  I  A W ) -DELK 

1 r  ( APRr < T HF ) -9o  )  1007,1007,1006 

T  HE  =  -  9  P 

GO  TO  1012 

F  I  Nn  A  ■T'-NNA  M  i '  6  F  h  F  OR  O^ER.  T  1  n  .\j  TP  RIGHT 
I  A  N|  T  r  K  *  J  A  W 

!  ■'  (AO  —  i  A  \!  T  )  1  009.1010,1010 

I  AMt  !  AJ  '  -  MCA 
T  aP  =  r  p,  *P*ANPA«  (  I  A  w  -  1  ) 

[-  (  a  R  c  f  t  T  M  E  )  -  R  ’  )  1012, 1012,1  l!ii 

T  J  E  =  ?  0 

C a |  r i •  i  t-  TEE  i;  ' E P i  t X  ANTE N A1  wp  T a r- 1-  O f  ThF  GIVFm  anTEmma 

c:j5tpe  =  o  r <  the /fcnead ) 

A  ■*  P  ;  'P''A'P»S!NT<'Ci'.  t-l*COST»-E  ) 

PHl-rn  'P?*rr)VTHp*wPH 

A  R  F  (  'A'  T  )  r  A  R  E (  I.  NT)*AEP*CCS(-mi  ) 

A  1  M  <  'A  1 T  ! s  A  I  P (  1  A  NT  )  *  A  M  P  •  S  I  K IFa  >  ) 

STOP*-  A'F  PeaSa  RIFP  TNFPRM..T|n„ 

P  h  I  u  r  I  \|  r  ,  |  w  A  '/  E  )  =  P--  I 
A  Mpu  f  I  '■  i  ,  [WAV*-  ):AMP 

C  3  N  T  I  M  m  f 
C  3  NT  '  AM,  P 

I  4  A  V  F  s  N  '  V  P  *  1 

I  ('■  W  V  -  -  1  W  A  \J  F  l  1  i  1  6 , 1  0  0  1  ,  1  :  O  1 

COAT  1  A  P 

H-«n  N;  T  >F  ANO  v C A f  N  1  N G  OEEP.  T|T,  INFORMATION 
R  -  a  n  q  1  1 ,  A  A  H  I  S  .  R  i\ 

R  r  A  r-  9"lP.^FrT,RCA  T,ICCAN,|«.AIJ,,RAM1 
CAI  r-i-i  T-  R  P  A  N  "  1  N  G  OPERATION  F 1 1  AJ  9  T  A  N  T  9 


mfnnMii— MUMinri.n  n -.i.iiMiiiniM.iiiiiiiia ait  — 1-iiniMt  r  ni  inm 


ooo  o  r:  a  ooo  odd  ooo  ooo 
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W  IDTh:'  ISAM-1  )  •  A  M  I 

SinPrSt-r  r-koIPTH/2 

I7  (  c  I  t  )  1 1M  .  i  151 , 1 1 S 1 

1150  S  I  nc  =  3(S(i  +  s  I  DF 

1151  CDNT!N"F 

D=l  T^SCanT^SAM  1/360 
A  \I0I  c  J=  .  7  n  7  *  A  N  C. ’ s 
C 

C  P  4  I  M  T  n  RFCORC  l  N  F  I  ’RM  A  T  I  ON 

c 

PRINT  0n?3.ANClS,RN,SECT,ISCAN,<?i'ANT,  lb  AM, SAMI 

c 

c  storf  h  record  information 
c 

A  \-A(14)=AN0IS 
A  -<eA  <  15  )  =RN 
4  A  ( 1  (S  >  =SFCT 
ARFA(17)=ISCAN 
A  RE  A (  1  5  )=SCAM 
ARFA<1<9)  =  ISAK 
A  RFA ( ?n  )  rSAM  I 

WRITF  H  RECORD  information  ON  to  tapf 
call  wPL^fPFAN, AREA, 50  ) 

B  *  C!  I  N  N  ’  N  G  OF  SCAN 
DO  1  ?5n  l  S  =  1 » I  Sr  AN 
STORE  s  RECORD  information 

A  RE  A  f  11=SH 
AREA  (  ?)=SFTMIM 
AREA!  3  )  =  I P  U  N1 
A  RF  A  (  4  )  r  N P(J N 
AREA!  S)=IS 

ZERO  max  OUTPUT  CATA  VARIABLES 

DATASM  =  0  .  n 
DATATMrO  •  0 

IMITIA  IZE  TIME 

T  =  <  I  S- 1 )  *SC  AN'T  -I'EL  T 

BEGINNING  OF  CAiCULATlONS  FOR  SAMPLE  POINT 

DO  1 ?2B  I  P  =  1 , ISAM 
C 

C  CALCULATE  Samplm  CENTER 

C 

C  r  N  7  r (  I  P - 1 ) • S A M I  4 S I DE 
I  7  f  rFVT-36fl )  3153,1152,115? 

115?  CFN'T  =  CCN  r-360 
1153  CONTINUE 
C 

C  ADVANCE  time 

C 

T  : T  +  r  Ei  T 


no~»  ono  ooo  noon  ooo  ooo  ooo  ooo  o  o  o  ooo  ooo 
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S<1P  W‘'V^  PHASE  STEP  CALCLl  A  i  I UNS  f  OH  MRS!  SCAN.  FIRSI  SAiPlE 
I  •*  (  r )  Il9n.1l9n.ii60 

H  -  0  t  N  N  1  N 1  >  OF  WA  ./g  PHASE  STEP  CALCULATIONS 

1160  03  117?  IWAVFsl .NWAVE 

I  '  WAV-  PMASF  STEP  sO.O  SKIP  CAI.i:UL  AT  IONS 
1*  (PSTPP{ 1WAVM  )  1161.1172.1161 
CAL  HCL '>TI-  PRESENT  H  AVE  PH  ASF 

1161  W3PH  =  Pc'Tl-P(  IRA  VP  )  •  1  /SCANT 

CALCM.ATH  NFW  A'TENNA  VOLT  ARK  f OP  AH  El  EHFNTS  IN  ARRAY 
03  1170  I  ANT « 1 , ' 0 A 

I*  (JMkh( IanT, 1  -  AVF  )  )  1169.1170.1169 

1169  P4lsPH'H(lANT,  i,AVP)*WPHH<  I  WAVE) 

P4?SPHH(IANT, 1aAVF)*WPPH 

APF<  I  A*  T  >  r  ARM  !  ANT)  ♦  AMPH  (  I  AM  i  .  !WiVb>«fL<'S<PA2)-C0S<PAl>) 
AjMMA'TJsATMlAM )*AMPH<  I  AMT  ,  I WAVF )*(S  N(bA2>-S!N(PAl>) 

1170  C3NT : N"F 

STOPF  PRt-SFMT  *iVE  PHASE 

1171  H’WHI  T  -  A'/F  )  sUPP* 

1172  CjNUN  'F 
1190  C'JNT  1  Ni'F 

C  A  L  P 1 1  ■  A  T  ir  MfirilLA  T  I  ='  s  FACTOR  FOR  i  I  MF  T 
a>i«;-(i  -^,:av)*pc  •p«‘'oS(ponk«t  > 

D p T C tM  |  m k  ANTCMI.A  Pl.FMPNT  TT  lHt  LEFT  TPF  S*MP  fe  CENTER  AND  THfc 
AMGi  f 9  RHTWFFN  'RF  SAMP|.F  CF’  tk R  AMT  AIMACFNT  ANTENNA  fclEHEMS 

K0=<fE  T-AWOA?) /ANPAM 

C3F|  cf''FNT-AN0A?)-«N0A»(KC-1  ) 

C  "•  F I  r  .  N  JA-rn^L 

CALC'  L  ■*. T *•  THE  I.  'NEAR  I NTgRPn.  A  1  !  (IN  FACTORS 

C3FI  a*i  n-|  /ANOA 
CO  El  AP  =  P  'Ft  P/AN  A 

I-*  KT  IS  NOT  fjfcThFKN  THF  LlMitS  1  -  NHA  .  CORRECT  THE  VALUE  OF  KC 

I-"  C<n  1194,1  194.1195 

1194  K:rNT  A ♦*(’ 

1195  I r  (KC-MOA)  1197.1197.1196 

1196  KCsKT-'RA 

L  f  F  T  P'Mi  F  PF  i I R  r  1  AND  RIRHT  RnNK  FOR  PHRs? 

1197  03  l?2n  iniosi.9 
C 

C  Z-RO  «''jv  ACCRV  LAIOP 
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C 

S  J*M  =0 
s jm  p  =  n 
c 

C  P^trnR«  JP-RATlON  f  OR  FACF  FlFRENT  IN  1hR  F  A  N  K 

C 

119*?  DD  1  P 1  7  iSF  =  1  ,NAS2 

G  3  rn  (i?nn,l?P4)  ,  I  D I R 
C 

C  F  !  MTI  £i  F'^FMS  I ' '  TnF  LFFT  AMI  FIGHT  [V  SCANNING  E '  t  M  E  N  I  Ibr 

C  FTP  TH-  FFT  SCANNING  RANK 

C 

120(1  1  -=KC-<  I  SF-j  ) 

C 

c  CAI  ri'LAT1-  T  HP  M  ■ '  I S  F  VOLTAGE  f  OR  THt-  AMH-NNA  ELEMENT  COMMON  To 

C  RDTH  T^F  I  FFT  Af'C  FIGHT  SCAMMN'G  ANTFNNa  BaMKS 

C 

I ■  (  I  S- - 1 )  1201.1201 ,1702 

1201  CAM  R  n(RN.RNl)R) 

AVRPRF  =  A'|0ISJ*HMCP 
CALI  R'vO(RN.RNDP) 

A'JIMRF  =  aNOISJ*R'iCP 
A  NRFPs  \NRFRF 

A  V 1 M P  s  a  N I MRF 
GO  T 0  1703 

c 

C  CALC'L’Tc  THF  Nn  IS-  VOLTAGE  *•  OR  THt  ANTFnna  ELEMENTS  TO  Tut 

C  LEFT  a  " n  RIGHT  OF  f HE  SCANNING  Fl.EMFNT 

C 

1202  AVRFPsAMREL 
A  V  I  MR  s  A  N  I  Ml. 

1  203  CALI.  R  O(Rim.RMjf) 

ANRFI  s  A  NO  I S J*RNOP 
Cali  R  « D  (  RN  ,  RN'DP  ) 

AMIM'  s^M  USJ*RN')P 
GO  TO  i?DR 
C 
C 

C  FIND  E.  FmFNTS  In  THE  LFFT  AML  RIGHT  OF  SCANNING  ELEMENT  ISb 

C  FOR  TH1-  RIGHT  So  A  NN  I  NR  BANK 

C 

1204  I.sKi+'Sr 

Hs  II  *1 
C 

C  CALCILAT-  THF  NniSF  VOLTAGE  *■  OR  THfc  ANTFNNA  ELEMENTS  TO  THt 

C  L-FT  A f  H  RIGHT  OF  THE  SCANNIN  G  Fi  EMFNT 

C 

I-*  (ISF-I)  1705.1205,1206 

1205  ANRFI  r  A  nRFRF 
A  M  I  mi  saNIMRF 
GD  T P  1707 

1206  ANRFI  sANHEP 
AMIM'  sANIMR 

1207  cam  rmtkrn.rndp) 

AVRFRsANO I S  J  *  R  N  n  P 
CALI  RNn(RN.RNDP) 

A  V  I M P s  A N O I SJ«RNPP 
C 

C  CORRECT  IR  AND  'L  I O  FIT  PETUEFN  THF  I  IMITS  1  -  NOA 

C 


ooo  ooo  o  o  o  o  ooo  ooo  ooo  non  ooo 
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120A  lr  (ID  l?ng,  1205, 1  210 
1200  I.sNDA*!1 

1210  Ir  (lL-N'Ui  1219, 1212, 1211 

1211  I _  *  I »  -NO* 

1212  Ir  <'R>  1213.1213. 1214 

1213  Hsn"A*TR 

1214  I-"  (IR-NO*)  1216.  1216.1215 

1215  M«I«?-M0A 

1216  C3NT  I  M > IP 

FlNn  ThF  VOl.  T  AGE  Of  THF  SCANNING  ELFMENT 

S9E*( A*S*AHF( ILMANRtl  >«CCEi  AP*(AHS#A»r(  |R)*ANRFR>*CDEl.A 
S  I  Ms  <  A  MR*  A I  M  < J  L )  ♦  A  f  '  l  HI. )*CCEl  A  P  ♦  ( A  M  S  •  A IM(|R)^aNJMR> *CU£L  A 

AJD  OM  THF  PHASE  OElAV  ANT  a  ,  •  !J  THE  VOLTAGE  TC  THE  BANK  ACCUMULATOR 

S  JM  f  sS"M  I  *SIM«CnSP  USE  )-SRE*S  IN  P  <  ISF  ) 

SJMR=SnMR*^  I  m«S  1  K  P  (  I SF  )  ♦SREsf  D^P(  I SE  ) 

1217  CONTINUE 

STORt-  THE  RANK  VCL'AGE 

S.IHRP(  I  0  i  R  )  sSUEH 
SJMJHt IDiPIsRUMI 

1220  C0NT I N"F 

CAL^HL-T*-  THF  VOLTAGE  SUMV  Of  TH*-  right  AN0  left  hank 

S  JMP  =  S"MPF(  1  >  ♦SnKRF-  <  2  » 

S  JMI  sS"H  I  M(  1  )*S'>H  [M(2  ) 

D»TaST=S.3RT  (SUHP««24.Mjm  f  *«2  » 

C  A  L  C  u  l  a  t  r-  THE  Vf'LTAGE  niFFEP-NTF  DETWEHN  THF  RIGHT  AND  LED  BANK 

SJMP  =  Shm-VF(P)-S'HRE(1  ) 

S  JM I  =Sl  H'M(2)-S"PIM(1) 

DATAnT=S'JRT(SU»*P«»2*SUMl»«2) 

CALCI'I  ate  THF  HHASH  DIFFEREni  F  RHTWFfcN  THF  RIGHT  AND  LEFT  BANK 
AMO  aO.'IIST  TP  FIT  pETwCFN  Tm-  LIMITS  -179.090  TO  IflO  OtbREt-S 

A  MGR  =  A  TANPCSHM  D  <  2  >  ,Sl)MRE<  2)  )*UONRAr 
A  V'Gl  sA  T4N2CCMM  I  -(  1  )  ,  StJMRfc  <  1  )  )  •( DURAT 
PATPrA  :GR-aNGL 

I-*  f 1  A  rp*i«o  >  1221  . 1221  .1229 

1221  OATPrO- TP*36P 
GO  T'l  1224 

122?  I-*  fl  A  IP-1  An  )  1224.1224,1223 

1223  OATPrpATH-360 

1224  continue 

STORE  THE  SCAf  OLThuT  data 

0  A  T  A  S  {  'PJsOATASl 
OATil  (  I o  )  a  P  A  T  A  D 1 
0  A  T  A  e ( IPIsPATP 
DATA'  <  'PJsCFET 

STORT-  S  WFCOPI  'NFiRMATION 


n  o  n  o  nnnn  non  n  n  n  n  nonnrsn  nnn  non 
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I  I  I  =  (  ! 

A  A (  T  I  T  )*CFNT 
A3Ca(  I  I  )«0A'  4ST 
A  3F  A  ( I  !  ! ♦? ) *PA  T  ACT 
A  3E  A  (  I  I I*J>*nAle 

FIND  T“>E  MA*lMU«  SUM  AND  DIFFERENCE  VALUE 

I"  (aR^FCDatasI  )-DATASM)  1226  .  1  226,  1  225 
1225  DAT  Ai-MiriATAST 

1228  Ir  (oacrtnATAm  )-DATAf!M)  1228 . 1  228. 1227 
1227  l)ATA'*M  =  nAT4nT 
122P  CDNTiN"c 

sropt  s  wprnRD  inf  irmation 
AREA!  *>=hatask 

A  3R  A  (  7)=n4TA(lM 

w a i t»-  <  xPnoRn  inf.rmtion  on  to  tape 

CALI  W">.*T<''FA\,Afifc«.«in) 
p3int  •>  oFcnRn  informatics 

P31NT  V019.SFTN(,M,  IRUN.NRLN 
P3!mT  0(128.  IS,  PaTa^M^aTALM 
P  3  I  Ml  2026 

P  3 1 N  T  S  RECORD  SAMH.E  PCIM  OUIPUT  TATA 
0  3 1  mT  8  ►  ILF  SAvPL‘  POINT  OlJ'-PUT  DATA 

DO  1240  iPsl.lSAM 

P3!MT  0n->7.nATA'-(If-),nATAS(IP),nATAP(TP),nATAP(TP), 

1  ( amack • f I ) . 1 1*1 .11 ) 

P^prnR-  iPF-R  A  T  IONS  F  OR  THF  I  i  *JP  PRIMpR  PLOT 

I l«( ( OAT-SM/I no*CA  I AS( IP)  )/0Al ASM) •50*1 
I  ?s<  <  ru  TM1M/1  8D*CA  l  Afl<  |P)  )/nATAn»)*50M 
I  3*  (  'PataP(  IP)*‘83.6)/^60)»8f'*1 
A.  I  went  t  *46 

P3  t  Mi  0-110  .  {  AI.  U  f  (  i  I  ) .  u  «1 ,81  ) 

A  .  |  ’  1  >»0 

A . In*  i  1 2  >*44 

bumt  OOIO.TAU'  EC  ’  I),  II*1.8i  ) 

A  I  MA  (  '  ?  i  s(* 

A_I\'M  1 >  * 24 

Pl|V»  2010,  TALI  F<  :  I  ),  I  1*1,8*  ) 

A  .  IM^-  f  1>sn 
1240  CiNT’M'-c 
1280  contim'ic 

6'0  fF  ruTA  RUN.  SiART  CALClII  A 1  ( ONS  OVER  Fr  R  NFXT  DATA  RUN, 

03  h  O’-'C  T APF  A f-'C  EXIT  PRCGRam  IF  LAST  r-ATA  RUN  MAS  BEEN  CuHflETEU 

I-’  r-R  in-IRUM  *995,8995.850 

T-*p  (I1-  IATA  RU'-S  "FF  THE  EM"  OF  TnF  TAPE,  PrUT  -RROR  STAlfehfeM 
Avn  tpwmimatr  t^ECi'Tlrw  of  pri-t-ram 
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a « o o  9  o  s  o 

g o  t i  «9R9 

c 

r  w9[tp  Hh-  T  F  fJ  M  1  < .  A I  i-lLF  linn  T  h  F  ENO  n^  T  H  I-  IAPF 

C 

8995  C  4 L  1  r  OSF ( DF  AN ) 

8996  CAU  n^F^(PFAN.nFAs, ITAPF.0.4HSJMU) 

A  9  F  \  (  i  )  -  F  h 

CALL  WRl  W  (OFAN  .  ARE  a  ,  ?  ) 

C 

C  C_n<;F  .MU  P  F  W  I  N  *  ThR  tape  AM'1  R  X  I  T  ThP  PROGRAM 

C 

8999  CALI  Ci  OSR  (  OFaN  ) 

CALI  SRTWN(  I  T  APR , 1  ) 

CALL  R  V  I  i 
C 

C  ••#*#*£#•#####*  £«*»•*##•*•#**#•****•***•****••******•**•*•••*••**• 

c 

C  T  -IF  r/n  FORMATS  ARP  AS  FOLLOWS 

C 

90  0  8  F  ^ r  < A  T  (315) 

9009  FORMAT  (  I  5  »  2F 1 0  .  2  »  15,  FlO.?) 

9010  FORMAT  (  l*5,3Flfi.«> 

9011  FORMAT  ( 5F 8 . 2  ) 

901?  FORMAT  (?Flf).2,?I5.F10.2) 

9013  FORMAT  (Fl 0  ,4,  i 5 ) 

9019  FORMAT  ( 1 R1 /5X , i 3HSFT  NUMRFR  =  .110.14H  •  FUM  HUMMER  ,l3»4h  OF  , 
IJI3.5W  R||MS) 

9020  FORMAT  (  /5X.25HANTRNNA  ARRA 

1Y  INFORMATION//  5>  ,  POHNLIFRFR  OF  ANTFNNAS  :  ,  5  7X  ,  F6 . 2/5X  ,  3  2P-D  I  S  T  A  NC 
2E  rpf'M  r-MTRB  Tu  SiRFFM  =,5x.r6.2,7M  MFTFRR/  5  X  ,  33HD  I  S  T  A  MCE  FROM 

3 S C R R >  M  TO  ANTFN'A  =,4x,F6.?.7H  MFTRRS//  5 X , 1 2R SC ANNER  MOUb// 
45X,?QMn||MRFR  OF  SC.MMMG  ANT-nNaS  =,BX,P6.?/  5X ,  16HC0PHASAL  ANGLb 
5=  ,  21 X . 1  6 . 2, RN  PFRGPFS//  5X,iAHWAVE  I  NF  I/C-M  A  T  I  ON  /  /  5X , 19HC APP 1 FR  ►  Fit 
6QjFwrv  =,1RX,F6.2.4P  MM7/  5X  .  22HMOIUIL  A  T  i  ON  FRRCUEniCY  = ,  1 5  X  ,  F  6 . 2 , 4h 

7  <H7/5X  .  23NPFRCI-N  T  OF  MOOLLAIICM  ■  •14X.-6.2/  5X.31HNUMBRP  OF  INC 

8  I  0  R  m  T  waVPFPONTv  s,6X,F6.?/> 

0021  FORM  *  T  (5X.4PRA''E.5X,4RWAVF,  <  (  6  X  .  4*W  A  VE  )  ,  6X  ♦  5HF  M  A  S  E /5  X  ,  4  HNUM  h  ,  5  X  , 
13HAMP,  7X.  5HPMASF.5X  ,3PAZM,  7x.  4HFi.EV.6X,  4  HSTEP/  ) 

9022  format  (SX, ?3,lv,5(4X,F6.2  )  ) 

9023  FORM'T  (  /<5  X  ,  1  7  H*  C  I  ^  E  I  mFOPM  A  I  I  CN//  5  X  ,  1  9HST  A  N  D  A  RD  DR  V  A  T  I  ON  »,16X,F 

16.?  /  5X.l3N9Nn  STARTER  S.24X.F6.2//  5  X  »  2  7H  AN  T  6  Nim  A  ARRAY 

?  SC  A  *  <-c  i  T!MPS//5X,l5PcFnOR  CFMTER  =  ,  ?2x ,  F  6 . 2 , 8H  DEGREES/ 

3  5X . i 8 -mumper  0<  9.ANS  =  ,19x,la.2/  5X  ,l5PS 

4C  A  N  Tier  ,?2X.F6.2,RP  9-C  /  5  X  ,  2  5HN  UM  HP  R  'F  SAMPLb  POINTS 

5r,15>,M.?/  5X.26WNCRI-MFAT  -PTw-EN  pOjMTS  =  ,  1 1  X  , 1  6 . 2 . 8h  DtGrtEES) 
9025  FORMAT  (  cX.'!3"''CAN  NUMPFu  =,13//  5x.14pMaX  SUm  DATA  s,2X,F6.2/ 

15X,1AM‘-ax  PfF  :  •  Ta  =,2y,Fp.?/) 

91126  FOOma.T  f /7X  ,  4PH-PF  ,  7>  .  XMRI  M,8X  .  3^0  |  f  ,  8X  ,  5RPHASF  /  7X  ,  5HS  1  GP  i  ,  6X  ,  4  H  U 
1  A  T  A , 2 t TV . 4  M  D  A  T  A  )  /  ) 

9027  FOPm-.T  (  2X  .  4  (3X  .F8.2)  ,5X,  AJ  .  A  ,  10  (4X,  A1  .4)  ) 

9030  FORMAT  (  1M+.MX  .51^1  .4) 

9050  FIOm-t  ( 1  Ml  /  5  X ,  3  7  H  T  H  P  [UTA  RL.NS  OFF  TME  FNC  CF  ThF  TAPb/ 

1  5  X  ,  pAHPPrGR'M  hi  I  L  l  NOT  RF  FXFCUTED 

C 

C 

E  -in 


^EDUCED  PROGRAM  LISTING  NUMBER  ONE 

THE  SASIC  OPERATION  OF  THE  PROGRAM 


VEC  »299.8 
PI  ■3.14159265 
TWOPI  ■  PI«2 
CONRAn«360/THOPI 


IRUN  *0 

READ  9008.NRUN, IlAPR.HODE 


914 


-GO  TO  (924, 914), MODE 
CONTINUE 


•  TAPE  MODE  2  -  Ol  D  TAPE 

•  ADVANCE  TAPE  TO  END  OF  OLD  DATA 

•  SETNIIH.SETNUH  or  LAST  DATA  SET 


- GO  TO  930 

l»*V24  CONTINUE 


♦  •  TAPE  MODE  1  -  NEW  TAPE 

♦  «  Initial I7E  tape 


SETNUM*0 
I  >*30  SETNUH«8EtnUH*1000 
V50  IRUN  «IRIIn*1-* - 


♦ 

♦ 

♦♦♦♦♦♦♦♦ 


♦ 

♦ 

*  *  ♦ 


i  n'lfilir 


READ  9009.NOA, A.D.NAS.CPH! 

REaD  90l0.NWAVE,EBFOC»rHEnM.PrAM 

N0A4  •NO A / 4 
ANOA  »36o/NOA 
ANOA?  aANOA/? 

NAS2  -NAS/? 
wl  •vEc/frSoc 
CONSTl"TWOP!#D/wL 
CONST?«TWnPI#A/WL 
CONW  ■TWnP!*r«fcQM*1000 


DO  1000  l ANT«l«NOA-* - 

AIM( I ANT)«0.0  1 

ARE ( I  ANT )  *0 . 0 

DO  1000  IWAVE«l,NWAVE-« - 1 

AMPHdANT,  IWAVE)«n.O  t 

loon  CONTINUE - 

I  WAVE  «1 

1001  CONTINUE-* - 

READ  9oil.WAMP.WPHI .A7M.ELEV.PRT  n 

WPH  »WPM I /CONRAO 

COSELV«COS(FLEV/CONRAD) 

COSAMP«CORELV»wAMp 
CONEl  «CONSTl«COSELV 
C0NE2  «C0MST2*C0SElV 


•  antenna  element  VOLTAGE  CAI  C'lLATIONS 


IwAVE  *IWAvF*1 

IF  (NWAVE-IWaVE)  1016.1001, 10m - *■ 

3016  CONTINUE 

READ  9gi3, ANOIS.Rm 
ANOISJ«.707*ANOIS 

READ  P012, SECT, SCANT, I  SCAN, I  SAM, SAM? 

WIDTH  ■( ISAM-1 > •SAm I 
SIDE  ■SECT-WIDTH/2 
IF  (SIDE)  1150,1161,1151 
1151  SIDE  •360+SinE 
1151  CONTINUE 

DELT  *SCANT»SAMl/^6n 


r 


NbDtlCFD  PROGRAM  LISTING  NUMBER  TWO 

phase  r>Ei ay  pai culationr 


CONCH l *COS(CpWl /COnRAD )»CnNST? 

CnSM  *CO«<  (  (MAS2-1  >»AMfiA*AMr)A?>/CPMRAr)> 

DO  990  ISE«l  .NAS2-* - 

PDFLAY"C0NCH!*(C0S<  (  ( j  SR -1 )  »Amo  A  ♦AMO  A? ) /OPmRaD  ) -CnSM ) 

cnsP( ISE )*cns<PDEi  a y > 

S!NP( ISE>«SIN(UELAY> 

999  CONTINUE - - 


(•EDUCED  program  listing  number  three 

ANTENNA  ELFMFNT  VOl  TAGE  CALCULATIONS 


K  •( A7h/AN0A)+1 

DELK  ■Azm-AnoA*<k-1) 

DELKP  «ANOA-DELK 

WPhh( iHaVF)*0.0 

PSTEP< IwAvE)«PST/C0NRAD 

DO  1014  lPIR«l,2-* - 

1002  DO  1«M3  lAU=l.N0A4-< - 

- GO  To  H003.1008),  idJR 

1003  IANT  *K*1-IAU 

- IT  (IANT)  1004, 1004. 100* 

1004  JAnT  "NOA^IAnT 

)005  THE  *ANOA*(l -I AW)-DElK 
- IE  ( ARSE ( THF ) - 90  )  1 007 , 1 007 , 1  0 0 A 

1006  THE  a-90 

1007  GO  To  101? 

100«  IANT  *K*IAw 

- IE  (NOA-IANT)  1009.1010,1010 

1009  IANT«IANT-N0A 

1010  THE  ■DElKP«-AnOA»MAW-1, 

- IE  ( ARSE ( THF ) -90  )  1012, 1012, 1M1 

1011  THE  *90 

1012  COSTHE*COS( THE /CONRAD) 

AMP  *COSAMP*SIN(CONE1*COSTHF) 
PHI  *CONE?*COSTHg*WPH 
APE (IANT) .ARE ( IANT) *AMp»COS( PHI ) 
AIm(IANT)«AJM<  JANT)*AMP*SlN(Puj ) 
PHIH( IAnT, IwAv6)*PHI 
amphuant,  I  WAVE  )  «AmP 

3013  CONTINUE - 

1014  CONTINUE— - 


i 


^EDUCED  PROGRAM  LISTING  NUMBER  FOUR 

TIME  DFPENOENT  PHASE  SHIFT  CAL  Till  AT  |  HNS 


- IF  I  T )  11P0.1190, 11*0 

1160  DO  1172  IwAvE*l,NwAvE —6 - 

- -IF  IPSTEPI  IWAvE) >  1161.1172.1161 

1161  WPPH  ■PSTEP( IHAVE)*T/SCANT 

DO  1170  UNt"T»MOA-^ - 

- IF  <AmPM<Unt.  IW*VP>  >  1169,1170,1160 

1169  PAl  *PHIH< I  ANT, IUAVE>*WPHH( IWAVE) 

PA2  «PHIH, I  ANT, 1WAVE)*WPPH 

ARE< I ANT>«ARE( I*NT)»*MpHl I  AMT, I WAVE > • < COS <P A? ) -Cqc ( p A 1 ) > 
A  I  Ml  IANT)«AIM<  IANTUAMPHU  ANT,  I  WAVE )  •  ( S  I  w  I P  A?  >  „S  I  N  <  PA1  )  ) 

1170  CONTINUE - i 

U7i  upmhiiMave)*wpph 

ii7?  continue- - 

1190  CONTINUE 


►'EDUCED  PROGRAM  LISTING  NIlMRpn  FIVE 

SCANNtNG  El  EMENT  AMR  SCANNING  OA*'K  VOLTAGE  ^Al  ClJALT  I  PNC 


KC  •  (CENT-ANOA?)/AnOAM 
CDEL  ,<CFNT-ANOA?>-ANnA*<KC-l  ) 

COELP  »ANDa -CDEL 
CDEL A  «CDFl  /ANOA 
cdclap«cdelp'anoa 

- IF  <KC)  H94,  U94,  ^95 

1194  KC  «N0A4KC 

1195  IF  ( KC "NO A )  1197,1197.1196 

1196  KC  ■KC-NnA 

1197  DO  1220  ID  IR«1 » 2  ^ - 

SUM  I  «0.0 

SUMP  *0.(1 

1199  DO  121 7  IsE*l,NAS?-<— - 

- GO  TO  <12no.l204), iniR 

1200  IL  ■KC-lsE*l 
IR  ■ I L ♦ 1 

♦  #  NOISE  SIGNAL  CaiCULaTION  F*r  LEFT  ♦ 

♦  ♦ 

<>0  TO  1208 
1204  IL  MCMSE 
IR  *IL*1 


♦  •  NOISE  SIGNAL  Cai CULaTIRNS  TOR  RIGHT  ♦ 

*  ♦ 

1208  IF  (II  11209, 1209,1210 

1209  IL  aNOA+lL 

1210  IF(IL-NOA)  1212.1212.1211 

1211  IL  aiL-NOA 

1212  IF  (IR)  121.1, 1213,1214 

1213  IR  «NO A* I R 

1214  IF  (IP-NOA)  1216,1216,1215 

1215  IR  aiR-NOA 

1216  CONTINUE 

SRE  »(AMS«ARE(II  >*ANREL)#CDFLAP4(AMS»APE(lR)4ANREP)#cnEl  A 
SIM  »(AMS»AIM(  Il»^ANIML)*CDFI  AP*(  6MS*A  I  MM  RUAN  J  MR)  #CDEl  A 
SUM  I  ■SUMI*SIM»C05P( ISE )-SPE»SlNP( ISE) 

SUMR  ■SUMR4SIMaSlNP(ISE)*SRE*COSP( ISE) 

1217  continue - » 

SUMREI IDIR)«SUmR 

SUM  IM(  1 D  I R )  »SI)M  I 

1220  CONTINUE  - 


57 

FnP  LF-rt  RAW 


^EDUCED  PROGRAM  LISTING  NIJMRf«  six 

determination  or  twf  noirf  signal  vm  take  - 


IF  < i sE-i >  1201, 1201.1202 

call  pND(RN.RMDP) 

ANRERFaANOIS.JaRNDp 

CALL  BNDIPN.RNDP) 

ANIMRF*ANO!SJ*RNDp 

ANRER  aANRERF 

AN  I  MR  ■ AN  J  MRF 

GO  T0  1203 

ANRER  «ANREl 

AN  I  MR  «ANI  Ml 

call  RNO(PN.RNDP) 

AnREL  "ANO I Sj#RnDP 
call  PNDIRN.RNUP  ) 

AN  I  ML  •ANniSj#RNDP 


REDUCED  PROGRAM  LISTING  NUMBER  SEV&M 

determination  of  Twe  notSf  signal  voltage  -  fdr  ritJh  ra“k 


- IF  (IRE-1)  1205.1205.1206 

1205  AnREL  *ANRERF 
animl  * an  I MRF 

- GO  TO  1207 

'*—1206  ANREL  «ANPER 

animl  *animR 

l-»“12o7  call  pND(pn.RNDP) 

ANRER  a ANO I S J#HNDP 
CALL  BND(RN,PNDP) 

ANIMR  * ANO I S j#RnDP 


DEDUCED  PROGRAM  LISTING  NUMBF*  F'®HT 
calculation  of  system  outputs 


SUmR  ■SUmRF(1)*SUmRE<2> 

SUMl  ■SUMIM(1)*SIIMIM(2) 
DATAST*SQRT  <  S1)HR»*2*SUM  !••?  > 


SUM*  ■SUmRF(?)*sUm*E<1 > 

Sum  I  "SUM I m ( 2 ) "Sum  I H ( l ) 

DATADT«SORT(SUMR»«!»*SUHl»*2> 


1221 


IS 

=122? 

1223 

1224 


angR  •aTan?(SUmim(2).s1)m«E(?i  >*CON«AD 
ANGL  ■ATAN?tSUMlM(l),SHMBEM  >  >*CONbAD 
n*TP  ■ ANGR- ANGL 
-IF  IDATP^ISO)  1221. 1221, 1??2 
DATP  *0 ATP*3AQ 
-GO  TO  1224  _ 

IF  (DATP-180 )  1224,1224,1223 

DATP  *OATp-3A0 
CONTINUE 


. 1  1  V 


mr 


APPENDIX  C 


SAMPLE  INPUTS,  FOR  GCAAS,  AND 
HARD  COPY  OUTPUT 
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S  4  MPl  F  INPUT  CONDITIONS  FCR  uCAAS 


COL. 


1  5 

l  n 

20 

30 

AO 

50 

4  4 

4 

4  4  4 

4  4. 

4  4 

* 

1 

1 

1 

12" 

1  AS  .  “511  6  .  on 

48 

2o.no 

? 

7 

.  o  o  o n  5.0000 

0  .  0  00  0 

I 

.00 

''.DO  90.00 

30.00 

n  .no 

.50 

10.00  85.00 

3^.00 

30.00 

0 

.  o  oon 

0 

90.00 

2.00  ? 

30 

2.00 

T-*F 

h a w n  COPY  OUTPUT 

FOR  THPSF  SAMPLE  INPUT 

IS  CONTAlNFr  IN  THE 

F 01  LOU  I  MO 

THREE  PAUFS 

OF 

60 
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SET  mhmhPiV  =  1  000  *  RUN  NUMRER  1  OF  1  RUNS 

ANTFNNA  ARRAY  INFORMATION 


NIMRFR  OP  ANTENNAS  = 

D I  STANCH  frpm  CbMTFF-'  TO  SCRE  r*'  = 

DISTA"CF-  FROM  scrfep  to  ANTENNA  = 


120.00 

145.50  mfTERS 
6.00  MHTFRS 


SCANNER  MOPE 


NUHPFR  or  SCANNING 

ANTENNAS 

: 

48.00 

C03HAPAI  angle  = 

20 .00 

DERGEES 

W  A  YF  INFORMATION 

CARRIFR  FREQUENCY 

s 

7.00 

MWZ 

MODUI  AT  1  ON  FRFQUFNCY  s 

5.00 

Kt-Z 

percent  of  modulation  = 

0  .00 

NIHRFF  nr  I  nc  I  pent 

WAVFFRONTS  a 

2.00 

WAYF  WAVF 

WAVE 

WAVF 

WAVF 

PHASE 

NURP  AMP 

PHASF 

aZm 

fcl  EV 

stfp 

1  1  .00 

0.00 

90 .00 

30.00 

0.00 

?  0.50 

10.00 

85.00 

35.00 

30.00 

NOISF  1 1  'PORMAT  I  ON 

ST  4  N  D  A  R  ti  DEVATIOM 

s 

0.00 

RND  STARTER  a 

0.00 

ANTENNA  ARRAY  SCAN 

SETTINGS 

SECTOR  OF'MTFR  = 

90 .00 

DFGFFES 

NIMRFP  Or  SCANS  = 

2.00 

SCAN  T  I MF  a 

2.00 

SEC 

NUMPFR  OF  SAMPLF  POINTS  a 

30.00 

INCRFMFNT  RFTWFEN 

POINTS  a 

2.00 

degrfes 

■-*■* - m 


SET  M I  i  M  h  p  =  1  000  *  RUN  NtMHER  1  OF  1  RUMS 

SC4M  MJMRt-R  =  1 


X  4 

4  4 

X  ♦  4  4 

4  ♦  N - MX - 4 - 

4  4  4 

X  4  4  4 

X  X 

- M  —  — - — - —  —  — - 

V  OVOO  fOonaOO  <0  M  'O  'O  N  V  N  H<OfM  H  O  ®  IA  fO  IT  N  fO  V  N 
.jj  4ro(\jNH'Oir»OMOroc\iin'Oo4HK5(\jinc>HcO'0'OOir\N 

c r>  <  . .  •  •  • 

<.  t-  oro^-r^oir\csjr^m»HC\jK)orvir»ooj'Oooorv5rfOr-iooo'C\iv 

i  <a  ocv^-'Or^ir»Of^Nmcvia‘OC\j  i  wno^isumo^^noojom^ 

Q.  Q  HrldHdHdHHHH  I  I  H  ddddddddd 

till  I  III  III 


KjlTi  D^CXOw'C^rf^  O  *  *  £  C  C  V  r'  r  C  X  ^  O  <1  n  r  r  X 

*occcr'o*c\)<vc«coof^*rva*i©^r*cvirr^oir*’<rcvr^ 


f^;in^;r^a)aDa'^dvr^rvr^r')oc\iajr'jrvfvmTHa)'0'or>‘'omrjc\jo 


~  < 
o  c 


c  O- 

a* 


o  fv 

rv  r- 


|V(^5r®^t--»f^5r'OLrOCD®'COrdK)P^fO^CVmfOOn-  <5  rH  C  H  H 

our^»HCJCvjvrNji^rvr;r^^doof';^dtdorv^oc\jf-jc>jvcucviGur)o*H 


u  u 

<i  < 


r  h- 
a)  c 


(\J(^(\iH^^'0<0'0NO0H4O4^f0'0f0if>Oinf0OH(\Jf0C^W 

*HrHC\JC\JCMC\lr4«d 


<  < 
c  c 


X  u. 
IT  C 


C.CC.COOCOGCOGCCOOOC.OCOOOC  C  O  O  C  CO 
OOOOOOOOOOOC  ODOCCOOODCDO  3  O  O  C  O  O 


X  X 
<  < 
r  x 


LL  a 
a  cr 
c  — 

rr  rn 


dWirNO'»-foir'Na  tHt-oiTKa^foirNr^aiHioirr^o* 
«c>'0'0>0'Cr,^r,^r^rsvrv  acoccoccccooooooooo'' 


ro  it  O' 

*d  tH  *d 
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N  AMk  STNPS) 

EQtiiP«CARDRE,PRlNTF  I 
SINGLE  (ARFA.OFAN)i 


PRO  t>  RAm  i  STHP«; 

STANDING  F  OR  |  SAMPi  E  TAPP  HANDLING  PROGRAM  SFGME-'T 
MR i t  Ten  RYi  ROY  p  hUNNINGmAUS 

WRITTEN  FORT  R AD  I OLOCA  T  I  UN  RESEARCH  LABORATORY 

oEPARTMFNT  OF  FLEOTRICAI.  FNGINPhR  I  Yf, 

university  of  Illinois 


INPIT  CaPDS  NEtESSaRY  FOR  OPFRaTION: 
one  PROGRAM  OPERATION  CARQl 

C 1 9 >  MSET  «  NIIMRFP  OF  oaTa  SETS  To  RE  LOCATED 

(15)  ITAPE  •  TAPE  DRIVE  NUMBER 

the  following  paros  are  needed  for  each  data  set  to  be  located 

DATA  SET  AND  RlfN  NllMPER  CARD  I 

< I 1 n )  SETNUM  •  Data  set  REFERENCE  number  to  Rt  LO0ATFD 

(15)  IRtlN  «  DATA  run  REFERENCE  NUMBER  To  BE  LOCATED 


PLALE  A|lL  ADDITIONAL  INPUT  CARDS  FOR  TNF  GIVFN  DATA  SET  -  PUN 
AT  TMIS  POINT 


DIMENSION  AREA(500l,DFAM650) 

6H«*H  F 
hH««H  H 
sm«ah  s 
DFAS«600 


READ  9000, NSET, ITAPE 

CALL  SFTBN( ITAPE, 1) 

CALL  FXBLK( ITAPE, D 
CALI  BLKLNUTAPE.O.O) 

T  SE  T  »  0 
GO  TO  101 

10"  CALL  SETPN( ITAPE, 1) 

101  iset.iset*i 

IF  (NSFT-ISET)  8500, 10?, 10? 

10'  PEAl  9oo1»SETNUM,  IPIIN 
ASSIGN  8000  TO  LAST 
GO  TO  10« 

101  CALI  CLOEE(DFAN) 

CALL  SETFD(DFAM) 

1QY  CALI  OPEMDFAN.DFAR, I T APE, 0 , 4MS T *U > 


CALI  RDlRIDFAN, AREA, l AST > 

10*  IF  < AREA<2)-SETNUM)  103,107.103 

10*  PRIM  9r>5o»SFTNUM 

107  NRUf  *AREA (4) 

IF  1  NRtJN"  I  RUN  )  lfl8.lll.lU 

10^  PRpT  R051,SeTmUm,  |PHm,nRUki 
10g  CALL  ClOSE(DFAM) 

(iO  TO  100 

no  CaU  RDLR(DFAN.  area. LAST  ) 

111  IF  <  IRUN-AREA<jM  110,112,110 

112  continue 

pRjnT  90<9,SeTnUh» IRUn 
jSCan«ArEA(17) 

!SAh  ■ AREA ( 19 ) 


r 

C  PLALE  ALL  OPERATIONS  FOR  THE  INPUT  INFORMATION  AT  THIS  POINT 

r 

r 

[’ 


DO  2999  lS*l, ISCAN 
CALL  RDLR1DFAN, AREA, l.AST) 


PLACE  ALL  OPERATIONS  FOR  THE  SCAN  AT  THIS  POINT 


2990  CONTINUE 


CALL  CLOSE(DFAN) 

00  TO  100 

8U0  CALL  CLOSE (DFAn ) 

CALL  SETPN( ITAPE.I) 
PRINT  9052 
CALL  EXIT 

85o  PRIM  9099 
CALI  EXIT 


900  FORNAT  (215) 

900“'  FORMAT  (110,15) 

9)4*  FORMAT  (lHl/5X.l?HeET  NUMBER  ,»i0/  *X.l2HRl|M  Ni  |MbF  p  .H.AH  FOu"0 
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1  ) 

90511  FORMAT  <lHl/5X,l?H9ET  NUMBER  ,110/  5X.17HN0T  FCIND  ON  T*PF  ) 

9051  FORMAT  (1H1/5X.12HSET  NuMBFR  ,110/  5X.1PRPUW  NU^PF0  (> 

1  4  if  H  <  GREATER  THAN  NUMBER  OF  Rll^S  I N  TWF  SFT  <,!3.1M)  ) 

9n5?  FORMAT  < //5x.2MTApE  OPERATION  0«  TAPE  F»ROR  ) 

909si  F  OR*  AT  (1H1/5X.3AHALI  OPERATIONS  COMPLETED  -  PPOiFA-  END  ) 

910n  FORMAT  <1H1/5X.1?HSET  NUMBER  #lJ0»14w  •  RUN  NUM-it-R  ,I3.4h  OF  . 

1  U,  5H  RUNS) 


end 


LISTING  OF  PROGRAM  STRO 


o  r. 
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PAGE 

001 

M 

name  stroi 

W 

EOll  1  PaCAHDRE 

,PRINTP| 

H 

c 

SINGLE  (AREA 

,  OF  An  >  1 

r. 

c 

progr  am  t 

stro 

r 

c. 

STANDING  for 

1  SIMULATOR  tape  read  out 

c 

r 

c 

WRITTEN  PYI 

POy  E  MuNNINGHAUS 

WRITTEN  F OR t 

RADIOLOCATION  research  laboratory 

C 

u 

r 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 

UNIVERSITY  OF  ILLINOIS 

L, 

COMPLETED! 

MAY,  1970 

r. 

r 

REVISED! 

oCTORER.  197fl  -  For  NEW  TAPE  HANDLING  PROCFDllRES 

* 

c 

0 

inplt  cards 

NECESSARY  FOR  OPERATION! 

i; 

one  program 

OPERATION  CAfiOi 

c 

0 

c 

<  1 5 >  iTape  *  Tape  dhiVf  n'Jmrer 

OlMfcNSlON  ARRAYS 

tJ 

OIMtNSjON  ALINE<51). AMA^htll) 

OlMkNSlON  AREAMlO>.nFAN<650> 

r 

r-  SET  OPERATION  CONSTANTS 

r 

EH*4H  e 
l)FAS*600 

t; 

REAL.  PROGRAM  OPERATION  INFORMATION 

r ; 

»EAL  9000* ITAPF 
r-  INITIATE  LINF  PRINTER 

00  900  11*1,51 

900  ALINECII)*0 

no  901  11*1,11 

901  AMAFK<!I>*27 

r 

C  OPEN  TAPE  ON  UNIT  NO.  ITAPF 

n 

y 

CALL  SFTRN( I  TAPE , 1 ) 

CALL  FYBLK(ITAPE.l) 

CALI  0LKIN ( I  TAPE ,0,0) 

ASSIGN  4000  TO  LAST 
GO  10  1001 
r 

r;  ADVANCE  TAPE  AND  PRINT  OUTPUT  DATA  UNTIL  TFRMINAL  FILE  IS  PEACHED 

1 0 o r  CALL  CLOSE(DFAN) 


i 
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PA(»E  0  0k 

CALI  SEtF  D  <  Of  Am ) 

1001  CALL  OpE  MPF  AN.  DF  AS, ITA  PE  *n,4nSjMU) 

CALL  RDlP(DFAn,ArEA.(  A*?T  > 

IF  UREA<1)-EH)  1100,3000,1100 

v  PRINT  H  RECORD  INFORMATION 

O 

1100  NRUN«AREA(4) 

(TO  7 O  1102 

1101  CALL  RDlR(DFAN, AREA. LAST) 

110?  NwAVE«AREA<13) 

IRI)N«AREA<3) 

PRINT  9019, ( ARFA( 1 1 ) » 

PRINT  9020, ( ARFAj I) ) . 11*5,13) 

PRINT  90?1 
kkk«21 

DO  H05  I  W  A  vE«l  ,  NW  A  vF 
KKK*»KKK*4 

PRINT  9022, I  WAVE, (AREA* 1 1 >, I J*KKK ,KKK4) 

1105  KKK*KKK*5 

PRINT  9n23,(ARE*Un,IH»lA.?0) 

I  SC *N« ARE  A ( 17 ) 

I SAM*APE A  <191 

r 

print  s  record  Information  for  each  scan  in  tme  data  run 

r; 

DO  1260  IS-l.ISCAN 

CALL  RDLPIOFAN, APEA.l  AST) 

1200  PRINT  9019,(ARFA(I|),||»2,4) 

DATASM*AREA<  6) 

DATADM«AREA(  7) 

PRjrT  9025, (ARFA( n )* I!B*»7) 

PRINT  9026 
DO  3250  IP.I.ISAm 
I  I  I  ■ « I  P-1 >  *4*11 
CENT  «aPEA<III  ) 

DATAST.AREA<  1 1 1*1 » 

DATADT«AREA(II  U2) 

DATAP  *AREA  < 1 1 1 ♦ 3 > 

PRINT  9027, CENT,DAT AST, DATADT, DAT AP, (AmARK( 1 1 ), I  I *1 , 1 1 > 

C 

C  PRENORM  OPERATIONS  FOR  LINE  PRINTER  PLOT 

11* < (DAT ASM/10 0+DAT AST) /DAT ASM )«50*1 
I2*( (DATADM/10flO*DATADT)/DATADM)*50^1 
13*  ( (DATAP*183.6)/3fS0)*50*l 
*L  I NE  <  1 1 )  «46 

PRINT  9030,(ALINE<II)*I!,1»51) 

ALINE( 1 1 » «0 
A  L I N  E « I2)«44 

PRINT  9030, (ALINE(II), II. 1,51) 

ALINEC 1 2 ) »0 
AL I N  E  < 1 3 ) «?4 

PRINT  9030, (ALJNE<II>, 11*1,51) 

ALINE< I3)*0 
125r  CONTINUE 

1260  continue 


°  age  ooj 
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GO  RACK  AND  PRINT  W  RECORD  I  NFORM  A  T 1  ON  FOR  THE  NCXT  DATA  Wl'N 
IF  LAST  DATA  RliN  IN  T  HE  DATA  SEt  MAS  RPEn  PRINTED,  AnVAMCf  T  APF 
To  NEXT  TATA  SfT 

if  »nriin-irun>  lnon, inno.iini 

TERMINAL  FILE  HAS  RFF N  REACHED#  PRINT  ALL  DATA  PRINTED  STATPHgMT 

3 1*  0  r  PRjnT  9070 
«0  TO  5000 

enl  of  tape  has  bfen  reached.  p*< i wT  error  statement 

4000  PR  I A  T  9060 
GO  T o  5000 

close  and  rewind  tape  and  exit  program 

5uor  call  Close <  df  an ) 

CALI  SFTPN(lTAPE.l) 
call  Exit 


The  I/O  FORMATS  ARE  AS  FOLLOWS 
9000  Fo^AT  (  i 5 ) 

911R  FORK  AT  (1H1/5X.13HSET  NUMBER  ■  ,110, 14H  •  RUN  NUMBER  ,  IJ.4H  OF  , 

1 1 3, 5h  R^S) 

9020  FORMAT  (  /t>X,?5HAMTENNA  ARPA 

lv  I'TORmaTion//  5X,20wnUMRFR  OF  ANTENNAS  «, 1 7X , F6 . 2/5X , 1?hD I  ST ANC 
2E  FROM  CFNTFR  TO  SCRFEM  ■  #  «5 X ,  F6  .  ?  .  7H  MFTERS/  5X ,  33«n  I  ST  ANCE  FROh 
3SCREEN  TO  ANTENNA  «,4X.F6.?.7H  MFTFRS//  5X, 1 2HSCANnFR  MODE// 

45y,  if  9HNIJMRER  OF  SCANNING  ANTENNAS  «,8x.F6.?/  5X ,  16HCOPH  AS  Al  ANGLE 
5*  #  ?1  X,  F6 . 2»  8H  OERGFES//  5X,16HWaVF  INFORMATION//  5  X  ,  19HC  A  RP  I  E  R  F^E 
6UUENCY  ■.1«X.FA.?,4H  MHZ/  5X . 22HM0DUL AT  I  ON  FREQUENCY  m. 15X . F6  .  ?, 4H 
7  KH7/5X.23HPERCENT  OF  MODULATION  * , 14 X , F6 . 7/  5X,J1HMUMREP  OF  PC 
8lDE*T  WAVEFRONTS  ■,6x,Fh.2/I 

9021  FORMAT  <5x#4hwAVE,5x.4HWAVF. 3 (6y , 4HW A VF > , 6y , 5HPH ArE /5y , 4HNUM0 , 5x , 
13hAmP#7X#5hPHASE»5X,3hAZm«7X,4hfLFV,6X,4hSTeP/) 

9022  FoPmAt  (5x» 1 3  #  1 X,5(4x#F6,2) ) 

9  n  2  3  FORMAT  < /5y ,  17HN0ISE  INFORMATION//  5x  ,  1  9MST  AND  ARD  DEvATION  «,lflX»F 
16.2  /  5X.13MPND  STaRTfF  »,24y,F6.2//  5 < , 27HANTENNA  ARRAy 

2  SCAN  SETTINGS//5X.15HSECT0R  CEnTEB  ■ . 9?X, F6 . 2, 8^  PEGrFES/ 

3  5X.18HNUMPER  OF  SCANS  ■  »19X,F6.2/  5X  ,15HS 

4CAN  TIME  ■  ,2?X,F(S.?,«H  sFC  /  5X ,  ?5HNLIM8EP  f  F  SAMPLE  POINTS 

5"»12X,F6.2/  5X.26HINCRFMENT  RETwFEN  POINTS  « ,  11X . Ef , ? , 8H  DEGREES* 

902’’  FOR’-AT  (  5X.13HSCAN  VUMPpR  ■ ,  f  3  /  /  5v,14mmAX  SJm  nATA  =,?X.F  P.?/ 
18X#  1 4HMAX  DIF  DATA  *,2*.F 8.?/) 

9ii2*  PQRfAT  (/7X.4HPnRE,7x,THSUH,«X,3HDTF,8Y,5HPHASE/  7x  ,  5HS I  GHt  ,  6X ,  4*1) 
lATA,2(7X,4HDATA>/> 

9U27  Format  f2x,4<3v#F8.?>,5x,A1.4,ln<4v,A1.4n 
9l|3r  FORMAT  (1H*,50X»51A1.4) 

9(i 6  Format  <1H1/5x.  36HtHP  END  of  the  TAPF  has  been  PEACHED) 

907*1  FOR*  AT  (1H1/5X,  38HALL  THE  DATA  ON  THE  TAPE  HAS  BEEN  READ) 


end 
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